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ABSTRACT 


We  set  out  to  design  an  extended  target  elassifieation  seheme  while  determining  the 
target’s  range-and-Doppler  loeation  with  the  use  of  adaptive  waveform  for  a  closed-loop 
cognitive  radar  platform.  To  that  end,  this  work  is  divided  into  three  objectives:  1)  in 
support  of  determining  range-Doppler  locations,  we  investigate  the  ambiguity  function  of 
the  matched  waveform  called  eigenwaveform,  2)  in  support  of  target  classification,  we 
look  at  an  adaptive  waveform  technique  called  probability-weighted  eigenwaveform 
(PWE)  and  introduce  two  new  waveforms,  and  3)  we  design  an  integrated  range-Doppler 
map  and  extended  target  classification  scheme. 

In  this  work,  we  show  that  the  fundamental  properties  of  ambiguity  function  for 
extended  targets  are  different  when  compared  to  classical  waveforms  for  point  targets. 
We  improve  on  the  adaptive  waveform  called  maximum  a  posteriori  PWE  and  introduce 
two  new  waveforms  called  match-filtered  PWE  and  two-stage  PWE.  We  propose  an 
integrated  range-Doppler  map  and  identification  scheme  for  multiple  moving  extended 
targets.  Performance  comparisons  in  terms  of  joint  probability  of  identification  and 
determining  targets’  range-Doppler  locations  with  traditional  wideband  waveform  and 
the  three  PWE-based  waveforms  are  shown.  It  is  shown  that  the  three  PWE-based 
waveforms  perform  better  than  the  classical  wideband  waveform. 
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EXECUTIVE  SUMMARY 


In  this  dissertation,  we  apply  the  optimum  waveform  design  for  a  deterministie  extended 
target  (also  known  as  eigenwaveform  [1])  to  ambiguity  function  analysis.  Interestingly, 
the  fundamental  properties  (peak  and  volume)  of  ambiguity  function  for  a  known 
extended  target  are  different  when  compared  to  ambiguity  functions  of  classical 
waveforms  for  a  point  target  [2],  Extended  target  responses  vary  greatly  from  target  to 
target,  from  one  aspect  angle  to  another,  etc.  For  the  purposes  of  illustration  in  this  work, 
we  form  two  distinct  target  types  (oscillatory  and  less  or  non-oscillatory)  and  analyze  the 
ambiguity  function  of  these  two  types  with  the  use  of  eigenwaveforms.  Ambiguity 
function  depends  also  if  the  transmit  waveform  is  one  pulse  as  in  continuous  wave  (CW) 
or  multiple  pulses  as  in  pulse-Doppler  radar  systems.  For  multiple  pulse  eigenwaveforms, 
increasing  pulse  repetition  frequency  results  in  increasing  Doppler  sine-lobes  separation 
but  increasing  pulse  repetition  frequency  results  in  the  reduction  of  the  unambiguous 
range  just  as  in  traditional  pulse-Doppler  systems.  However,  due  to  the  increased  time 
interval  of  the  return  (with  the  use  of  eigenwaveform)  Doppler  resolution  improves  at  the 
expense  of  range  resolution. 

Based  on  optimum  waveform  design  for  a  deterministic  extended  target,  we 
investigate  and  expand  the  notion  of  the  adaptive  waveform  called  probability  weighted 
eigenwaveform  (PWE)  that  is  used  for  target  recognition  with  a  cognitive  radar 
platform  [3].  PWE  is  designed  and  implemented  for  non-moving  target  identification 
problem  in  a  cognitive  fashion  where  the  waveform  weights  of  eigenwaveforms 
corresponding  to  an  ensemble  of  targets  to  form  the  transmit  waveform  is  updated  via 
previous  return  (received  measurement)  and  likelihood  calculations  under  the  constraint 
of  transmit  energy  or  number  of  transmissions.  We  propose  two  new  PWE-based 
adaptive  waveforms  in  addition  to  an  existing  one.  The  three  adaptive  waveforms  are: 
maximum  a  posteriori  PWE  (MAP -PWE),  match-filtered  PWE  (MF-PWE)  and  two-stage 
PWE  (TS-PWE),  which  adaptively  uses  likelihood  information  from  previous  and  present 
stage  to  form  next  waveform  in  a  specific  way.  It  is  shown  that  the  three  PWE-based 
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waveforms  perform  better  than  the  classieal  wideband  waveform.  It  is  also  shown  that 
MF-PWE  and  TS-PWE  seem  to  perform  better  than  MAP-PWE. 

Our  ultimate  goal  in  this  dissertation  is  to  design  an  integrated  algorithm  that 
addresses  the  combined  problem  of  multiple  moving  target  identification  while 
determining  accurate  range  and  Doppler  locations.  Combining  range-Doppler  map 
(RDM)  technique  with  the  PWE-based  adaptive  waveform  techniques,  we  propose  an 
integrated  range-Doppler  map  and  target  identification  scheme  for  multiple  moving 
extended  targets.  The  overall  probability  of  correctly  locating  an  extended  target  in 
range-and-Doppler  and  identifying  target  type  is  clearly  a  function  of  received  signal-to- 
noise  ratio,  number  of  transmissions  L ,  PWE  schemes  and  the  target  responses 
themselves  via  maximum  target  eigenvalues.  In  practice,  a  priori  assumption  may  come 
from  other  sensors  or  some  intelligence  information.  At  the  end  of  this  work,  various 
performance  comparisons  in  terms  of  overall  probability  of  identification  and  probability 
of  correctly  determining  range-Doppler  locations  for  the  three  PWE-based  waveforms  are 
evaluated  against  a  classical  wideband  pulse  waveform.  All  the  waveforms  are  used  in  a 
closed-loop  radar  system.  It  is  noted  the  MAP-PWE,  ME-PWE,  and  TS-PWE  perform 
better  than  the  wideband  waveform. 
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I.  INTRODUCTION 


There  are  three  main  related  objeetives  in  this  work.  The  first  is  to  investigate  the 
application  of  matched  waveform  design  (also  called  eigenwaveform)  [1]  to  ambiguity 
function.  In  other  words,  we  are  interested  in  the  target  range  and  Doppler  implications  of 
the  eigenwaveform.  The  second  objective  is  to  extend  the  adaptive  waveform  called  the 
probability  weighted  eigenwaveform  (PWE)  first  proposed  in  [2]  and  investigate  its 
classification  performance.  As  such,  we  are  also  interested  in  identification  of  extended 
targets.  And  finally,  the  third  objective  (which  utilizes  the  results  of  the  first  two)  is  to 
design  an  integrated  scheme  for  the  combined  problem  of  range -Doppler 
location/detection  with  extended  target  type  identification  with  the  use  of  a  closed-loop 
radar  system  employing  adaptive  waveforms.  Such  a  system  is  an  example  of  a  cognitive 
radar  (CR). 

In  this  work,  we  utilize  adaptive  waveforms  based  on  eigenwaveform  which  is 
also  known  as  SNR-based  illumination  waveform  design  [1].  Thus,  there  is  a  need  to 
evaluate  ambiguity  functions  of  radar  systems  employing  eigenwaveforms  and  their 
properties.  We  investigated  eigenwaveform-based  ambiguity  functions  with  the  use  of 
three  types  of  targets  (oscillatory,  non-oscillatory,  and  one  based  on  an  actual  Boeing 
aircraft).  We  found  that  the  AF  properties  (peak  and  volume)  for  eigenwaveforms  to  be 
much  larger  than  that  of  traditional  waveforms  for  point  targets.  Unlike  traditional 
transmit  waveforms  whose  responses  totally  dictate  the  shape  of  the  ambiguity  function, 
both  matched  illumination  waveform  and  extended  target  response  contribute  to  the 
shape  of  the  ambiguity  function.  In  other  words,  range  and  Doppler  resolutions  are  not 
just  functions  of  the  transmit  waveform  but  of  the  target  response  itself,  which  makes  for 
interesting  ambiguity  functions.  Various  range-Doppler  trade-offs  are  made  with  the  use 
of  pulsed  eigenwaveforms. 

For  a  traditional  radar  system  where  targets  of  interest  are  very  far  in  range,  a 
good  and  common  model  is  to  assume  that  the  targets  are  point  targets.  In  additive  white 
Gaussian  channel,  the  received  waveform  is  therefore  the  scaled  transmit  waveform  plus 

noise.  The  probability  of  detection  is  then  a  function  of  the  received  waveform  energy 
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regardless  of  its  shape  [2],  Moreover,  the  ambiguity  funetion  is  dictated  by  the  shape  of 
the  transmit  waveform.  Waveform  designs  in  consideration  of  the  ambiguity  function  for 
point  target  model  has  a  rich  literature  [3]-[5]  and  these  waveforms  have  well  known 
properties  in  terms  of  range  resolution,  Doppler  resolution,  and  probability  of  detection. 
Our  interest  however  is  extended  targets,  i.e.,  targets  that  have  certain  impulse  responses 
(i.e.,  have  finite  time  support)  and  therefore  the  returns  do  not  just  depend  on  the  transmit 
waveform  but  rather  depend  on  both  the  transmit  signal  and  the  target’s  response  via  their 
convolution.  Given  constraint  in  energy,  it  is  shown  in  [6]  that  the  eigenwaveform  is  the 
waveform  that  maximizes  the  signal-to-noise  ratio  in  additive  white  Gaussian  or  colored 
noise.  In  signal-dependent  interference  however,  the  energy  spectral  density  (BSD)  of  the 
optimal  waveform  is  derived  in  [1].  In  [I]  and  [6],  the  emphasis  is  to  derive  optimum 
waveforms  for  extended  targets  while  in  [7]-[9],  the  emphasis  is  to  apply  the  adaptive 
waveforms  in  a  target  recognition  or  classification  using  cognitive  radar  (CR)  [10].  The 
traditional  targets  in  radar  field  are  assumed  to  be  point  targets  and  have  a  rich  literature 
of  waveform  design  and  detection  theory  in  [11]-[14].  However,  for  an  extended  target 
of  our  interests,  the  range  resolution,  ambiguity  function  and  probability  of  detection  are 
clearly  affected  by  the  target’s  impulse  response.  We  investigate  and  evaluate  the  range 
resolution,  ambiguity  function,  and  probability  of  detection  of  radar  that  employs 
eigenwaveform. 

Moreover,  we  evaluate  the  target  identification  performance  of  PWE-based 
waveforms  in  target  identification  problems.  To  this  end,  we  propose  improvements  to 
the  maximum  a  posteriori  probability  weighted  eigenwaveform  (MAP-PWE)  adaptive 
waveform  design  used  in  target  recognition  with  a  cognitive  radar  platform.  The  two 
improvements  are  called  match-filtered  PWE  (ME-PWE)  and  two-stage  PWE  (TS-PWE). 

Various  target  detection  and  classification  schemes  are  reported  in  many  articles 
and  a  few  are  referenced  here  for  the  interested  reader  [15]-[21].  But  our  research  is  to 
investigate  the  properties  of  ambiguity  function  with  extended  targets  and  develop  an 
integrated  detection  and  classification  schemes  in  a  cognitive  manner  [1],  [2],  and  [7]. 
Some  other  works  use  likelihood  or  correlation  characteristic  [22]-[26].  However,  our 
method  is  to  look  at  the  detection  probability  improvement  with  PWE-based 
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transmissions.  As  had  been  shown  before  [1],  [7],  [8],  and  [18],  the  detection  and 
classification  performance  of  the  adaptive  waveforms  are  better  than  just  employing  a 
traditional  wideband  pulse  waveform. 

Recall  that  our  ultimate  goal  in  this  dissertation  is  to  design  an  integrated 
algorithm  that  addresses  the  combined  problem  of  multiple  moving  target  identification 
while  determining  accurate  range  and  Doppler  locations.  Combining  range-Doppler  map 
(RDM)  technique  with  the  PWE-based  adaptive  waveform  techniques,  we  propose  an 
integrated  range-Doppler  location  and  target  identification  scheme  for  multiple  moving 
extended  targets.  Various  performance  comparisons  in  terms  of  overall  probability  of 
identification  and  probability  of  correctly  determining  range-Doppler  locations  for  the 
three  PWE-based  waveforms  are  evaluated  against  a  classical  wideband  pulse  waveform. 
All  the  waveforms  are  used  in  a  closed-loop  radar  system.  It  is  noted  the  MAP-PWE, 
ME-PWE,  and  TS-PWE  perform  better  than  the  wideband  waveform. 
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II.  RANGE  RESOLUTION  FOR  EXTENDED  TARGET 
ILLUMINATED  BY  AN  EIGENWAVEFORM 


For  a  point  target  illuminated  by  a  basie  pulse  (one  that  is  not  eompressed),  the 
range  resolution  is  easily  derived.  Since  our  target  of  interest  is  extended  (one  that  has 
finite  extent)  then  the  range  resolution  provided  by  a  basic  pulse  including  the 
eigenwaveform  has  to  be  analyzed. 

A,  RANGE  RESOLUTION 

For  point  targets  with  the  use  of  basic  one-pulse  waveforms,  the  minimum  (range 
or  time  delay)  separation  of  two  targets  is  dictated  by  the  length  of  transmit  waveform.  It 
is  defined  as  the  minimum  time  separation  needed  so  that  the  return  pulses  do  not  overlap 
each  other,  i.e.,  it  is  the  minimum  corresponding  range  separation  required  such  that  two 
point  targets  can  be  resolved  in  range.  This  is  easily  illustrated  in  Figl.a. 


Transmit  waveform 

ts 

n  t 

‘y' 

?2 

Return  waveform 

4-^ 

4-^ 

Return  waveform 


Return  tgl 


Return  tg2 


< -  2T  - ► 
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Figure  1 .  Target  range  separation  setting;  (a)  target  range  separation  for  point  targets, 
and  (b)  target  range  separation  for  extended  targets. 
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Let  L  be  the  time  duration  of  the  transmit  waveform.  For  the  return  waveforms 


not  to  overlap  eaeh  other,  two  point  targets  need  to  be  separated  by  7"^  =  L  /  2  as  shown 

in  Fig.  la.  For  basie  pulsed  waveform  (not  ehirp  or  waveforms  employing  eompression), 
the  eorresponding  minimum  range  in  whieh  two  target  ean  be  resolved  is  given  by 

R,=cT^=cTJ2,  (1) 

where  c  is  the  speed  of  light  (in  free  spaee).  For  basie  one-pulsed  waveforms,  minimum 
range  separation  is  usually  referred  to  as  range  resolution. 

In  praetiee,  range  resolution  ean  be  improved  by  waveform  design  (using 
eompression  teehniques)  that  ean  potentially  inerease  the  bandwidth  ( 5 )  of  the 
waveform.  As  sueh,  range  resolution  is  usually  given  by 

R^=cT^=cl2B.  (2) 

Using  minimum  separation  requirement,  if  the  target  has  a  response  with  duration  T  and 
that  the  eigenwaveform  has  same  time  duration  =T),  the  falling  edge  of  one  target 

return  has  to  be  separated  by  F  /  2  from  the  leading  edge  of  another  target  return  as 
illustrated  in  Fig.  lb.  Thus,  the  minimum  time  needed  in  order  for  two  extended  targets  to 
be  resolved  is  given  by  =  1  •57’ .  Thus,  the  minimum  range  separation  sueh  that  two 
targets  ean  be  resolved  with  the  use  of  the  eigenwaveform  is  given  by 

R^=cT^=\.5cT.  (3) 

For  an  extended  target,  both  the  eigenwaveform  and  the  target  response  dietate 
the  effeetive  bandwidth  of  the  return  signal.  As  sueh  a  simple  equation  eorresponding  to 
(2)  in  terms  of  bandwidth  does  not  easily  apply.  This  is  also  the  big  differenee  between 
point  target  and  extended  target  where  waveform  bandwidth  is  easily  manipulated  if  a 
target  is  truly  a  point  target.  It  is  elear  that  for  traditional  waveforms  (with  large 
bandwidths)  the  resolution  promised  by  (2)  will  not  be  realized  if  used  for  extended 
targets  (espeeially  if  a  target  has  a  narrowband  frequeney  response).  This  is  beeause  the 
frequeney  response  of  the  return  signal  is  the  result  of  multiplieation  of  the  frequeney 

responses  of  the  transmit  waveform  and  target  response.  Sinee  a  praetieal  extended  target 
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doesn’t  have  an  idealized  flat  frequeney  response  like  a  point  target,  the  frequency 
response  of  the  return  signal  will  surely  be  different  from  the  frequency  response  of  the 
transmit  waveform.  As  such  large  bandwidths  may  be  decreased. 

For  traditional  waveforms  with  time  duration  with  the  point  target  assumption, 
the  minimum  Doppler  separation  such  that  two  targets  of  different  speeds  can  be  resolved 
is  usually  taken  to  be  the  reciprocal  of  the  time  duration,  i.e.,  A/’^  -  HT^.  If  we  use  the 
same  definition  for  an  extended  target  with  finite  extent  T ,  it  follows  that  the  Doppler 
resolution  is  ls.f^=\l {T^  +  T) .  For  an  eigenwaveform,  it  is  therefore  A/'^=l/2r.  In 

other  words,  while  the  eigenwaveform  range  resolution  suffers  due  to  the  increased  echo 
return,  its  Doppler  resolution  improves  due  to  the  length  of  the  return,  which  accentuates 
the  time-frequency  resolution  trade-off. 

B,  REVIEW  OF  EIGENWAVEFORM 

In  Bell’s  work  [6],  he  derived  the  optimum  transmit  waveform  matched  to  an 
extended  target  in  continuous-time,  which  is  later  named  eigenwaveform.  For  transmit 
signal  x{t)  and  target  response  h{t) ,  let  s{t)  be  the  convolution  of  transmit  signal  and 
target  response.  Assume  a  matched  filter  in  the  receiver  for  signal  detection.  The  output 
of  matched  filter  is 


X{t)  =  j  s{t)s*{t  -  (4) 

-00 

where  r  is  the  time-delay  (location)  and  is  the  Doppler  shift  of  a  certain  return.  The 
s{t)  is  the  convolution  result  of  transmit  signal  and  target  response. 

For  now,  we  may  assume  zero-Doppler  since  the  energy  in  (4)  is  not  affected  by 
the  Doppler  shift.  We  may  assume  unit  energy  for  the  extended  target  in  our  signal 
models.  Any  other  energy  value  (point  target  or  otherwise)  can  be  accommodated  by 
scaling  the  energy  of  target  response  .  The  eigenfunction  {t)  corresponding  to  the 
maximum  eigenvalue  is  used  as  the  transmit  waveform,  aka  eigenwaveform  [1]. 


7 


Therefore 


40  =  ^max(0,  (5) 

and 

00 

^max^max  (0  =  J  ^„ax  (0^(P  "  r)d pcl T,  (6) 

-00 

where  /?  is  a  “dummy”  variable  and  L{t)  is  the  autoeorrelation  function  given  by 

CO 

L{T)=\\H{ffe^^^^^df.  (7) 

-00 

By  substituting  (6)  and  (7)  into  (4),  which  is  the  energy  of  convolution  of  transmit 
signal  and  target  response  becomes 

00 

=  j  ^max(^)^max^max*(^')«^^' 

(8) 

=  X  \  q  (T)q  {T^dz  =  X  E 

max  J  7/ max  V  max  max  ’ 

-00 

where  E^  is  clearly  the  transmit  energy  of  the  eigenwaveform.  That  is,  E=E^  . 

Compared  to  traditional  point  target  whose  E^  is  merely  E^  due  to  convolution  with  delta 
function  regardless  any  target  waveform,  E^  for  an  eigenwaveform  illuminating  an 
extended  target  is  amplified  by  X^^ .  Maximum  eigenvalues  from  various  target  types  are 
different,  which  means  that  the  energy  received  is  different  from  target  to  target. 

To  gain  an  insight  into  eigenwaveforms,  we  consider  two  complex-valued 
extended  targets  illustrated  in  Fig.  2,  which  are  arbitrarily  generated  for  the  purposes  of 
illustration.  In  the  top  panel,  the  magnitudes  of  the  target  frequency  responses 
(normalized)  are  shown.  Both  targets  are  of  unit  energy.  The  magnitudes  of  the  frequency 
responses  of  eigenwaveform  are  shown  in  the  middle  panel.  The  bottom  panel  shows  the 
magnitudes  of  the  frequency  responses  of  the  return  signals. 
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Target  1  Frequent^  Response 


Digital  Frequency 
Frequency  Response  Magnitude  of  Transmit  Signal  1 


Target  2  Frequency  Response 


Frequency  Response  Magnitude  of  Transmit  Signal  2 


Frequency  Response  of  Return  Signal  2 


Digital  Frequency 


Figure  2.  Top  panel:  Target  magnitudes  in  frequeney  response.  Middle  panel: 
Frequency  response  magnitudes  of  transmit  signal.  Bottom  panel:  Frequency 
response  magnitudes  of  return  signal. 


It  is  clear  from  Fig.  2  that,  the  eigenwaveform  distributes  most  of  its  energy  to  the 
target’s  most  resonant  frequency  band.  We  may  think  of  this  as  the  transmit  waveform 
capturing  the  most  of  the  energy  of  the  target  signal  via  the  frequency  domain. 
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III.  PROPERTIES  OF  THE  AMBIGUITY  FUNCTION  FOR 
EXTENDED  TARGETS  AND  THE  EIGENWAVEFORM 


In  this  section,  we  apply  extended  target  and  eigenwaveform  method  to  the  study 
of  ambiguity  function  (AF).  In  classical  ambiguity  function  studies,  the  received  signal  is 
assumed  to  be  a  replica  of  the  transmit  waveform.  In  other  words,  targets  are  assumed  to 
be  point  targets.  Again,  AF  analysis  of  various  popular  waveforms  such  as  pulsed,  phase- 
coded,  and  chirp  waveforms  assuming  a  point  target  is  mature.  Aside  from  the  excellent 
texts  cited  above,  we  point  out  a  few  select  ones  such  as  [27]-[31].  Our  interest  however 
is  extended  targets  and  how  various  transmit  waveforms  affect  ambiguity  functions.  More 
specifically,  our  interest  is  how  both  an  extended  target  and  its  corresponding 
eigenwaveform  design  affect  the  ambiguity  function. 

In  radar  application,  the  ambiguity  function  features  the  output  of  the  matched 
filter  as  a  function  of  Doppler  shift  and  time  delay  and  is  given  by 


I  x{t)x*  (t  - 


(9) 


where  x{t)  is  the  transmit  waveform,  is  the  Doppler  shift,  and  r  is  the  time  delay. 

Notice  that  this  definition  where  squared-magnitude  of  the  matched  filter  response  is  used 
is  from  [5].  Sometimes  the  magnitude  of  the  matched  filter  response  is  used  as  in  [4].  To 
avoid  confusion,  we  will  choose  one  and  use  the  squared-magnitude  definition  of  the 
ambiguity  function.  Again,  it  is  clear  from  (9)  that  the  received  signal  is  the  exact  replica 
of  the  transmit  waveform  (assuming  no  noise  or  interference)  where  the  target  is  assumed 
to  be  a  point  target.  Unfortunately  (9)  does  not  work  when  the  target  has  a  finite  extent. 
In  the  case  of  an  extended  target,  the  ambiguity  function  is  not  only  a  function  of  transmit 
waveform  but  is  also  a  function  of  the  target  response.  If  the  target  response  is  given  by 
h{t)  then  the  ambiguity  function  is  given  by 
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(10) 


\x{^Jdt 

00 

=  I  {x{t)*h{t)'\x[x{t-T)^h{t-Ty[ 

-00 

»  2 

-00 

where  s(t)  is  the  eonvolution  return  of  transmit  signal  x(t)  and  target  response  h(t)  . 

It  is  elear  in  (10)  that  for  extended  targets  both  transmit  waveform  and  target 
responses  dictate  the  ambiguity  function.  Our  interest  is  to  form  the  ambiguity  function  of 
the  eigenwaveform  given  an  extended  target  response  and  compare  that  function  with 
AFs  produced  by  traditional  waveform  such  as  wideband  pulse  and  rectangular  pulse 
waveforms. 


A,  DISCRETE  OR  VECTOR  SIGNAL  MODELING 


Due  to  the  advent  of  arbitrary  waveform  generators,  modem  waveforms  are 
designed  in  discrete-time.  Moreover,  modern  waveform  design  are  usually  simulated  by 
modern  computing  that  utilizes  discrete-frequency  techniques.  Proper  digital-to-analog 
conversion  (DAC)  easily  converts  the  waveform  into  practical  continuous-time 
waveform.  Assuming  proper  sampling  rate  through  Nyquist  theory,  we  can  utilize 
discrete-time  and  discrete-frequency  models  to  form  a  discrete-time  discrete-frequency 
version  of  the  ambiguity  function  for  extended  target,  which  is  given  by 


XV^n'Jdnl  = 


A^-1 


^5[z]x/[z-rjxlf, 


fdJ 


(=0 


I; (x[z]  * /z[z]) X  (x[z  -  r  J  * h\i  -  r„ ])*  x 
/=0 


(11) 


where  is  the  discrete  delay  and  is  the  Doppler  spread,  x ,  h  and  s  are  the  discrete 
transmit,  target  response,  and  return  signals.  N  is  the  length  of  return  vector  s  ,  =  0, 

1,  ...,  N-\,  and 
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In  general  the  length  of  may  be  desired  to  be  different  from  N .  That  extension  is 
easily  incorporated  into  (11). 


The  peak  of  the  ambiguity  function  is  located  at  zero-delay  and  zero-Doppler 


(/^«  =0)-  That  is, 


^5[z]X5*[z] 


X(xn*M/])x(x[zr/zm)* 


(13) 


It  turns  out  that  this  peak  is  not  constant  given  an  extended  target  considering  various 
transmit  waveforms  (in  which  the  proof  is  discussed  later).  Before  going  into  the 
extended  target  case,  let  us  examine  the  classical  point  target  case.  If  the  target  is  a  (non¬ 
fluctuating)  point  target  then  h  has  an  impulse  response  (i.e.,  no  extent)  and  thus  the 
peak  is  given  by 


X\.\  '■>  fdn  ] 


peak,  point 


=  E^  =  El 

s  h 


^x[z]xx[z]’ 


i=\ 


=  ElEl, 


(14) 


where  E,^  is  the  energy  associated  with  the  point  target  and  it  is  clear  from  (14)  that  the 

AF  peak  is  always  a  constant  no  matter  what  kind  of  waveform  is  used.  In  fact  for  point 
targets  (non-fluctuating),  the  volume  under  the  ambiguity  function  for  all  possible 
transmit  waveforms  is  also  constant  and  is  given  by  EIeI  .  If  the  received  point  target 

energy  and  transmit  energy  are  normalized  [E]  =1',EI  =1)  then  the  peak  is  1.  The 
properties  that  AF  peak  and  volume  are  constant  regardless  of  transmit  waveform  are  the 
basis  of  designing  waveforms  on  how  to  manipulate  ambiguity  sidelobes  in  AF.  We  will 
show  in  this  work  that  for  an  extended  target,  various  AFs  arise  depending  on  the 
transmit  waveform  and  the  point-target  properties  do  not  hold  for  extended  targets. 
Moreover,  we  discuss  how  we  accurately  implemented  our  own  ambiguity  function  (in 

MATLAB)  and  how  we  verified  our  results  in  the  Appendix. 

13 


B,  PROPERTIES  OF  THE  AMBIGUITY  FUNCTION 

1 .  Property  1 :  Given  an  extended  target,  the  peak  of  the  ambiguity  funetion  is 
different  for  different  transmit  waveforms. 

In  other  words,  the  AF  peak  is  not  eonstant  for  all  waveforms  i.e.,  in  (13)  is 

not  eonstant  unlike  in  the  point  target  ease  whose  peak  is  eonstant  as  given  by  (14).  We 
will  now  prove  Property  1 . 

For  eompaet  representation,  let  h  be  the  eomplex-valued  diserete-time  target 
veetor  of  length  n  ,  then  we  ean  easily  form  h  =  where  is  the  energy  of  the 

target  response  and  the  normalized  target  response  h  elearly  has  an  energy  of  1 .  Let  x  be 
an  arbitrary  transmit  signal  veetor  of  lengths  .  From  [7],  we  ean  form  a  normalized  target 
eonvolution  matrix  given  by 

h  0  0  0 
.  h  0  0 

h 

where  H  is  a  (2n-l)  by  n  matrix.  Thus,  any  arbitrary  transmit  waveform  x  yields  a 
eorresponding  ambiguity  funetion  whose  peak  is  given  by 

x{k-,l.]E  =  wo  •*[/])• 

/=! 

=  (x*h)''x(x*h)|'  =|(x*7^h)"x(x*7^h)|'  (16) 

=  El  |(Hx)''  X  (Hx)|"  =  El  |x"H"Hx|"  . 

It  is  elear  in  (16)  that  various  peaks  are  possible  depending  on  what  x  is.  In  faet,  it 
should  be  elear  that  the  value  of  the  peak  depends  on  the  how  x  eorrelates  with  the  target 

response  (represented  by  If  various  peaks  are  possible,  then  it  is  our  goal  to  find 

the  ambiguity  funetion  (from  all  possible  x  )  whieh  yields  the  maximum  peak  among  all 
possible  peaks. 
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The  largest  peak  can  be  found  by  maximizing  the  AF  peak  for  all  possible 
transmit  waveforms  x  given  H  via 


where  we  let 


R,  =  H^H. 

Based  on  the  structure  of  H ,  Ra  will  have  the  form 


i 

z  ■ 

i 

■  z 

i 

i 

i 

I*, 

i  / 


(17) 


(18) 


(19) 


where  all  the  diagonal  components  have  the  same  energy  as  h  (which  is  1).  It  can  be 
shown  with  some  effort  that  the  matrix  Ra  is  conjugate  symmetric,  Toeplitz  and  positive 
definite.  That  is, 


—  - H -  - H -  „  — H 

R,  =  H  H  =  (H  H)"  =  Ra 


(20) 


The  Ra  matches  exactly  the  definition  of  sample  correlation  function  from  [32]  but 
without  the  normalization.  We  may  refer  to  Ra  as  the  normalized  target  response 
autocorrelation  matrix. 


Notice  that  the  only  design  parameter  in  (16)  is  the  transmit  waveform  x  .  In  other 
words,  maximizing  the  peak  in  (16)  is  equivalent  to  maximizing  the  argument  x^Rax, 
which  is  given  by 


max£’? 

VxlH  * 


(21) 


15 


Notice  the  argument  inside  the  squared  magnitude  in  (21)  is  an  eigenvalue  problem  [11] 
where  the  eigenvalue  and  eigenveetor  relation  is  given  by 


and 


R/,q,  =  A.q. 


(22) 


A,.  = 


qfq, 


=  qfR/,q,, 


(23) 


where  Z.  is  any  eigenvalue  of  the  autoeorrelation  matrix  R*  and  q,  is  the  eorresponding 
unit-energy  eigenveetor  of  length  n  . 

Thus,  the  largest  peak  of  the  ambiguity  function  is  achieved  by  taking  the 
eigenveetor  q,^^  corresponding  to  the  largest  eigenvalue  as  the  transmit  waveform. 
In  other  words,  the  AF  that  produces  the  largest  peak  is  the  return  that  convolves  the 
eigenwaveform  with  the  target  response.  Sinee  q^^^  is  unit-energy,  we  need  to 
incorporate  any  transmit  energy  constraint  E^,  i.e.,  the  transmit  waveform  is 

X  =  ^/^q^ax  •  Thus,  the  maximum  AF  peak  given  a  target  response  eorrelation  matrix  is 


=maxA; 


Vx\H 


=ElEl 


qLR'-q 


peak  y^|jj 
2 


x^Rax 


=  /l"  EIe\ 

max  h  X  ’ 


(24) 


whieh  says  that  the  maximum  peak  is  given  by  the  squared-product  of  maximum 
eigenvalue,  target  energy  and  transmit  energy.  Moreover,  this  maximum  peak  comes 
from  the  ambiguity  function  produced  by  using  the  eigenwaveform. 

In  conclusion,  given  any  specific  target  response,  the  maximum  peak  value  of  the 
ambiguity  function  is  attained  by  taking  the  eigenwaveform  corresponding  to  maximum 
eigenvalue  as  transmit  signal  instead  of  any  arbitrary  waveform  (and  this  includes 
traditional  pulsed  or  wideband  waveform).  The  maximum  peak  value  is  Z^^^E^E^  as 

opposed  to  E^  in  (14)  whose  value  depends  on  both  how  well  the  transmit  waveform 

matehes  with  the  target.  Notice  that  a  peak  of  E^E^  can  be  attained  by  a  special 
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waveform.  This  AF  peak  is  attained  by  an  idealized  impulse  waveform  [7].  In  our 
research,  we  implement  the  wideband  waveform  as  practical  implementation  of  the 
impulse  waveform. 

2.  Property  2:  Given  an  extended  target,  the  volume  of  the  ambiguity 
function  is  different  for  different  transmit  waveforms. 


In  other  words,  the  AF  volume  is  not  constant  for  all  transmit  waveforms.  We  will 
now  prove  Property  2.  The  volume  of  the  extended  target  ambiguity  function  is  given  by 


F 


AF 


N 


fdn 


(25) 


It  can  be  shown  by  Parsavel’s  theorem  for  discrete  Fourier  transform  that 


/=0 


(26) 


Thus,  the  AF  volume  becomes 

1  N-\  N-\  , 

/=0 

N-\  N-\  - 

(A^— 1)  /=0 


(27) 


By  substituting  ij  for  i  and  substituting  4  for  ( ^  “  G  )  where  4  and  4  cover  the  range 
from  0  to  ( A  - 1 ),  the  total  volume  of  the  ambiguity  function  is  given  by 


N-\  N-\ 


f:..=ZZh[6]x/[4]| 

i2  =0  4  =0 

N-\  N-\  ^ 

=(i;kol  )x(i;k[4]  )=ei 


i,=o 


G  =0 


=  El 


x^R/iX 


(28) 


where  we  use  (22)  and  (24)  complete  the  last  line  of  the  proof.  It  is  clear  that  the  AF 
volume  for  extended  target  is  not  constant  and  various  waveforms  can  produce  various 
volumes  depending  on  the  extended  target  response  as  well  as  the  transmit  waveform. 
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One  of  the  waveforms  we’ll  use  in  our  work  is  the  wideband  (impulse)  waveform. 
An  extended  target  illuminated  by  the  wideband  waveform  has  a  volume  that  is  given  by 


=  At 

AF  n 


A/. 


dn 


N 


'L'LW.'jJ  =EiEi 

fdn 


(29) 


This  is  very  interesting  since  this  is  the  same  volume  produced  if  we  illuminate  a  point 
target  with  any  transmit  waveform.  In  other  words,  for  extended  targets,  only  the 
wideband  waveform  can  produce  this  volume.  With  the  use  of  the  eigenwaveform,  the 
largest  AF  volume  possible  is  given  by 


V,,  =  Ar„ 


A/. 


dn 


N 


fdn 


(30) 


In  other  words,  the  volume  is  the  same  value  as  the  largest  peak.  This  volume  is 
amplified  by  compared  to  that  of  the  volume  produced  by  the  wideband  waveform. 
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IV.  AMBIGUITY  FUNCTIONS  OF  ONE-PULSE  WAVEFORMS 


A,  EIGENWAVEFORM,  WIDEBAND  WAVEFORM,  AND  RECTANGULAR 

WAVEFORM 

Traditionally  for  point  targets,  a  basic  pulse  (usually  shaped)  is  a  common  choice 
for  transmit  signal.  For  extended  targets  in  this  work,  we  will  choose  two  pulse-type 
waveforms  and  compare  that  to  the  eigenwaveform.  The  waveforms  are:  wideband 
waveform  (i.e.,  a  very  narrow  pulse  compared  to  the  length  of  the  extended  target)  and 
rectangular  waveform  (i.e.,  a  rectangular  pulse  with  duration  less  or  equal  to  that  of  the 
extended  target).  Given  a  target  response,  we  derive  and  illustrate  ambiguity  functions  for 
these  two  waveforms  and  compare  them  to  the  ambiguity  function  of  the  eigenwaveform. 
In  Fig.  3,  we  illustrate  the  finite-length  of  the  three  waveforms  (wideband,  rectangular, 
and  eigenwaveform)  as  they  interact  with  an  extended  target  whose  response  has  finite 
support  T . 


Wideband  Waveform  Extended  target 

Convolution  >l« 


Return  signal 


I 


A 


A 


<-  T  ^ 

T  —> 

<r  T  —> 

OD 

Rectangular  Waveform  Extended  target 

Return  signal 

r  ^ 

1  1  Convolution  ^ 

A 

<-  T  -> 

<-  T  —> 

<-  T+r  — >• 

(b) 


Eigenwaveform 


Extended  target 


hd 


Convolution  ^ 


jA 


(c) 


Return  signal 


Figure  3.  One-pulse  transmit  waveform  for  extended  target:  (a)  wideband  waveform,  (b) 

rectangular  waveform,  and  (c)  eigenwaveform. 
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Let  be  the  time  duration  of  any  transmit  waveform.  As  mentioned  earlier,  we 
assume  the  extent  of  the  wideband  pulse  is  mueh  less  than  the  target  response  duration. 
As  sueh,  in  Fig.  2a  we  may  idealize  the  wideband  pulse  to  be  an  impulse  ~  0).  For 
reetangular  waveform,  the  pulse  width  is  T^-r  as  shown  in  Fig.  2b.  The  length  r  may 
be  thought  of  as  time-on-target  (ToT)  and  the  ratio  of  to  T  is  defined  as  the  time-on- 
target  ratio  (ToTR)  and  is  given  by 

ToTR  =  ^.  (31) 

In  Fig.  2e,  the  eigenwaveform  is  shown  to  have  the  same  length  as  the  target 
duration  and  thus  the  ToTR  is  equal  to  1.  For  a  wideband  waveform,  ToTR  =  IT 

where  T^«T  .  In  the  example  simulations  (in  this  work)  using  the  wideband  waveform, 
T^IT  is  1/31.  For  this  seetion,  we  eoneentrate  on  one-pulse  waveforms  to  gain  insight 

into  the  ambiguity  funetions  of  these  three  waveforms  as  they  eonvolve  with  target 
responses.  For  illustration  in  this  work,  we  form  two  targets  with  very  different  frequeney 
responses  and  set  both  the  transmit  energy  and  the  target  energy  to  be  unit-energy  {E^  =  \ 
and  =1).  We  will  keep  using  these  two  targets  for  eomparison  to  illustrate  the  faet  that 

unlike  in  the  ease  of  point  targets,  the  target  responses  play  important  roles  in  the 
formation  of  ambiguity  funetions.  In  Fig.  4,  we  present  the  magnitudes  of  frequeney 
responses  of  two  different  targets  (top  panel),  the  eorresponding  frequeney  responses  of 
the  eigenwaveforms  (middle  panel),  and  the  frequeney  responses  of  their  return  eehoes 
(bottom  panel).  From  Fig.  4,  we  ean  eonelude  that  ehoosing  the  eigenveetor  with  the 
maximum  eigenvalue  of  target  autoeorrelation  matrix  as  the  transmit  signal  is  tantamount 
to  ehoosing  a  band  of  frequeneies  (Fig.  4  middle  panel)  where  the  eehoes  (Fig.  4  bottom 
panel)  guarantee  the  largest  returns  in  terms  of  magnitudes  [1],  [7].  Interestingly,  the 
magnitudes  of  the  return  frequeney  responses  exhibit  suppression  of  the  less  dominant 
frequeney  bands  of  the  target  response.  This  turns  out  to  be  important  sinee  this 
suppression  effeet  will  translate  to  sidelobe  suppression  in  the  ambiguity  funetion  of 
target  returns  illuminated  by  eigenwaveforms.  For  the  purposes  of  illustration  in  this 
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work,  we  form  two  target  types  where  eaeh  represents  something  practically  meaningful. 
Target  1  (a  name  that  we  will  now  use  for  consistency  throughout  the  work)  illustrated  in 
Fig.  4  (left  top  panel)  with  two  dominant  bands  represents  targets  with  oscillatory 
tendencies  in  time  domain  while  Target  2  (also  to  be  used  consistently  herein)  in  Fig.  4 
(right  top  panel)  with  a  low-pass  shape  represents  targets  exhibiting  less  or  non- 
oscillatory  tendencies  in  the  time  domain.  Target  2  is  actually  a  low-pass  shaped  pulse  in 
time  domain,  which  represents  a  set  of  target  responses  that  have  a  strong  initial  return 
but  decays  off  in  time.  In  other  words,  we  have  two  target  responses  that  represent 
potential  practical  target  responses  and  thus  we  can  illustrate  how  each  extended  target 
type  will  respond  to  the  three  waveforms  we  are  to  investigate. 


Target  1  frequency  response  Target  2  frequency  response 


Normalized  Frequency  Normalized  Frequency 

Figure  4.  Top  panel:  target  frequency  response,  middle  panel:  eigenwaveform 
frequency  response,  and  bottom  panel:  return  signal  frequency  response. 

Moreover,  the  actual  radar  range  profile  of  a  Boeing  737-500  [33]  is  converted  to 
a  practical  target  response  and  used  an  example  to  generate  a  practical  eigenwaveform 
and  corresponding  AF.  This  response  is  shown  in  Fig.  5. 
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Figure  5.  Example  of  a  range  profile  of  a  Boeing  737-500  used  as  Target  3.  The  radar 

range  profile  is  at  the  bottom.  From  [33]. 


B.  EIGENWAVEFORM  VS.  WIDEBAND  AMBIGUITY  FUNCTIONS 

A  quiek  and  very  insightful  trade-off  eomparison  is  to  use  an  eigenwaveform  and 
the  wideband  waveform  as  transmit  waveforms.  In  this  seenario,  we  utilize  Target  1  ( n  = 
31  samples).  In  Fig.  6,  the  ambiguity  funetions  (3D  and  eontour  plots)  of  wideband 
transmit  waveform  and  eigenwaveform  are  caleulated  and  illustrated. 
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Target  1  Ambiguity  Function  (Wideband  waveform) 


Target  1  Ambiguity  Function  (Eigenwaveform) 
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Figure  6.  Ambiguity  function  and  AF  contour  plots  of  Target  1  illuminated  by 
wideband  waveform  vs  eigenwaveform.  Top  panel;  3-D  view,  and  bottom  panel; 

contour  plots. 


It  is  clear  that  the  ambiguity  function  of  eigenwaveform  is  more  compact  in  terms 
of  Doppler  (frequency)  response  compared  the  wideband  waveform.  The  cleaner  Doppler 
response  is  the  result  of  two  factors.  The  major  factor  is  the  fact  that  eigenwaveform 
exhibits  narrowband  characteristics  that  suppresses  the  less  dominant  target  frequencies. 
The  minor  factor  is  the  fact  the  length  of  the  return  echo  for  the  eigenwaveform  is  about 
twice  that  of  the  wideband  waveform.  Increasing  length  of  transmission  results  in  the 
increase  of  Doppler  resolution.  While  Doppler  resolution  was  defined  earlier  in  terms  of 
time  of  return  echo,  a  related  concept  defines  Doppler  resolution  to  be  half  the  effective 
mainlobe  width  in  the  Doppler  domain  given  a  fixed  delay  (usually  when  delay  is  zero). 
This  is  also  true  in  range  resolution  where  delay  (or  range)  resolution  is  usually  taken  to 
be  half  the  effective  mainlobe  width  in  the  time-delay  domain  given  a  fixed  Doppler 
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frequency  (usually  =  0).  Sidelobe  suppression  in  the  ambiguity  function  is  usually  a 
desired  design  characteristic  for  traditional  waveforms  for  point  targets.  Here,  the 
eigenwaveform  has  an  inherent  feature  of  being  able  to  suppress  sidelobes  in  the  Doppler 
domain  of  the  ambiguity  function.  Since  the  length  of  the  return  for  the  eigenwaveform  is 
approximately  twice  of  that  of  the  wideband  waveform  (recall  Fig.  3  and  Fig.  6  bottom 
panel)  the  mainlobe  length  (when  looking  at  the  delay  spread)  of  the  eigenwaveform  AF 
is  larger  compared  to  the  mainlobe  length  of  the  wideband  AF.  In  other  words,  the  range 
resolution  of  the  wideband  waveform  is  slightly  better  than  the  eigenwaveform. 
However,  the  eigenwaveform  AF  contour  is  clearly  much  cleaner  (due  to  suppression 
effect  and  the  Doppler  resolution  effect  of  longer  return)  than  the  wideband  waveform. 
Concentrated  energy  in  the  origin  of  the  AF  is  desired  in  radar  waveform  design  for 
avoiding  ambiguities.  Also  bigger  sidelobes  can  cause  a  threshold  to  be  crossed,  which 
can  cause  increase  in  probability  of  false  alarms  (PFA).  Compared  to  the  wideband 
waveform,  the  reduction  in  range  resolution  is  a  small  price  to  pay  considering  better 
Doppler  resolution  and  the  qualitative  sidelobe  reduction.  Indeed,  it  is  a  small  price  to 
pay  when  we  consider  the  peak  value  of  eigenwaveform  AF  compared  to  wideband  AF 
(since  peak  value  translates  to  improved  detection  performance).  In  Fig.  6  where  Target  I 
is  used  =  I ),  the  peak  of  the  eigenwaveform  AF  is  33.52  ( =  5.79  )  while 

the  peak  of  the  wideband  AF  is  1 . 

C.  RECTANGULAR  WAVEFORM  AF  FOR  OSCILLATORY  TARGET 

(TARGET  1) 

The  ambiguity  function  from  the  use  of  rectangular  transmit  waveform  depends 

on  r  relative  to  T  (i.e.,  ToTR)  in  Fig.  3.  In  this  section,  we  consider  two  rectangular 

waveforms:  a  low-ToTR  and  a  high-ToTR.  Using  a  high-ToTR  of  0.64,  we  illustrate  in 

Fig.  7:  a.  AF,  b.  AF  contour,  c.  Target  1  frequency  response  (magnitude),  d.  transmit 

signal  frequency  response  (magnitude),  and  e)  return  signal  frequency  response 

(magnitude).  Using  a  low-ToTR  of  0.09,  we  illustrate  in  Fig.  8:  a.  AF,  b.  AF  contour,  c. 

Target  1  frequency  response  (magnitude),  d.  transmit  signal  frequency  response 

(magnitude),  and  e.  return  signal  frequency  response  (magnitude).  The  AF  peak  of  high- 

ToTR  in  Fig  7  is  lower  (0.0298)  compared  to  the  AF  peak  (0.2482)  of  the  low-ToTR 
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rectangular  waveform.  This  is  because  transmit  spectrum’s  main  lobe  (sine-function  in 
frequency  domain)  shown  in  Fig.  7d  barely  overlaps  the  two  dominant  frequency  bands 
of  Target  1  (in  Fig.  7c)  thereby  reducing  the  energy  return  as  shown  in  Fig.  7e.  Notice 
that  there  are  also  some  subpeaks  and  sidelobes  in  AF  for  the  high-ToTR  rectangular 
pulse.  On  the  other  hand,  low  ToTR  has  a  better  peak  value  in  the  ambiguity  function  but 
a  qualitatively  worse  Doppler  response.  This  is  because  the  low-ToTR  rectangular 
waveform  has  a  wider  response  (in  frequency  domain)  than  the  high-ToTR  rectangular 
waveform  and  therefore  allows  for  more  frequencies  to  appear.  To  conclude,  given  unit- 
transmit  energy  =  and  unit-target  energy  =  1 ),  we  note  that  the  AF  peak  values 

for  both  rectangular  waveforms  are  both  lower  than  1.  Recall  that  the  wideband 
waveform  AF  has  a  peak  value  of  1  while  the  eigenwaveform  has  a  peak  value  of  , 
which  is  33.52  for  Target  1.  In  other  words,  rectangular  waveforms  (low-pass  shape  in 
frequency  domain)  are  not  good  waveforms  for  Target  1  when  it  comes  to  AF  peaks. 


(a)  Ambiguity  Function 

(Rectangular-transmit  waveform.  High  ToTR=0.64) 


0.03 


-100  -0.5 

Delays  in  samples  Doppler  in  Normalized  Frequency 

(b)  Contour  plot 


(c)  Target  1  frequency  response 


(d)  Transmit  signal  frequency  response  (Rectangular) 


Figure  7.  Ambiguity  function  analysis  of  rectangular  waveform  with  Target  1  (ToTR  = 
0.64):  (a)  ambiguity  function,  (b)  AF  contour  plot,  (c)  Target  1  frequency 
response  (magnitude),  (d)  transmit  signal  frequency  response  (magnitude),  (e) 
return  signal  frequency  response  (magnitude). 


25 


(a)  Ambiguity  Function  (c)  Target  1  frequency  response 

(Rectangular-transmlt  waveform.  LowToTR=0,09)  4i  i  >  i  '  i  i  ' 


0#l«ys  in  samples  Normalised  Frequency 

Figure  8.  Ambiguity  function  analysis  of  rectangular  waveform  with  Target  1  (ToTR  = 
0.09):  (a)  ambiguity  function,  (b)  AF  contour  plot,  (c)  Target  1  frequency 
response  (magnitude),  (d)  transmit  signal  frequency  response  (magnitude),  (e) 
return  signal  frequency  response  (magnitude). 


D.  RECTANGULAR  WAVEFORM  AF  FOR  NON-OSCILLATORY  TARGET 

(TARGET  2) 

In  this  section,  we  consider  a  less-oscillatory  target  or  non-oscillatory  target  in 
which  Target  2  is  a  good  example.  Again  we  let  E^=\  and  =1.  From  previous 
example,  we  have  already  gained  the  insight  that  the  wideband  AF  peak  is  1  and  that  it 
has  a  slightly  better  range  resolution  than  the  eigenwaveform  and  rectangular  waveform. 
Thus,  we  are  able  to  make  comparison  with  rectangular  waveform  AF  without  generating 
the  wideband  AF.  While  the  rectangular  waveforms  (low  and  high  ToTR)  do  not  work 
well  for  oscillatory  targets,  the  rectangular  waveforms  actually  perform  well  in  terms  of 
AF  peaks  for  the  non-oscillatory  Target  2  compared  the  wideband  waveform.  For  this 
example  scenario,  it  is  instructive  to  just  use  one  of  the  rectangular  waveforms  to 
compare  with  the  eigenwaveform.  We  choose  to  look  at  the  high-ToTR  rectangular 
waveform  AF  and  compare  that  to  the  eigenwaveform  AF  for  Target  2.  Using  a  high- 
ToTR  of  0.64,  we  illustrate  in  Fig.  9:  a.  AF,  b.  AF  contour,  c.  Target  2  frequency 
response  (magnitude),  d.  transmit  signal  frequency  response  (magnitude),  and  e.  return 
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signal  frequency  response  (magnitude).  Using  the  eigenwaveform,  we  illustrate  in  Fig. 
10;  a.  AF,  b.  AF  contour,  c.  Target  2  frequency  response  (magnitude)  of  target  response, 
d)  transmit  signal  frequency  response  (magnitude)  of  rectangular  waveform,  and  e.  return 
signal  frequency  response  (magnitude).  Notice  that  the  high-ToTR  rectangular  waveform 
AF  has  a  decent  peak  (232.38).  In  hindsight,  such  a  result  may  not  be  surprising  since 
Target  2  is  a  low-pass  shaped  response  and  that  a  rectangular  pulse  should  return  a  large 
echo  specially  if  the  durations  of  target  and  rectangular  waveform  are  comparable. 
Looking  at  rectangular  waveform  AF  contour  in  Fig.  9b,  the  overall  sidelobe  level  is  also 
low  (of  course  this  is  very  much  influenced  by  the  fact  that  both  target  and  transmit 
waveform  response  are  concentrated  near  the  zero  Flerz  frequency).  Notice  however  that 
AF  peak  resulting  from  the  eigenwaveform  is  still  higher  (377.45,  which  is  the  square  of 
/Imax  =19.4281).  Looking  at  the  eigenwaveform  AF  contour  in  Fig.  10b,  the  Doppler 

spread  is  narrower  than  that  of  Fig.  9b.  In  this  case  of  non-oscillatory  target,  the 
rectangular  AF  peak  is  clearly  much  higher  than  the  wideband  waveform  AF  peak  (which 
is  1). 


(a)  Ambiguity  Function 

(Rectangular-transmit  waveform,  High  ToTR=0.64) 


(c)  Target  2  frequency  response 


(d)  Transmit  signal  frequency  response  (Rectangular) 


Figure  9.  Ambiguity  function  analysis  of  rectangular  transmit  waveform  with  Target  2 
(AF  peak  value  =  232.38):  (a)  ambiguity  function,  (b)  AF  contour  plot,  (c)  Target 
2  frequency  response  (magnitude),  (d)  transmit  signal  frequency  response 
(magnitude),  and  (e)  return  signal  frequency  response  (magnitude). 
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(a)  Ambiguity  Function  (Eigen-transmit  waveform) 


iOO-, 


Delays  in  samples  Doppler  in  Normalized  Frequency 


(b)  Contour  plot 


(c)  Target  2  frequency  response 


(d)  Transmit  signal  frequency  response  (Eigenwaveform) 


Figure  10.  Ambiguity  function  analysis  of  eigenwaveform  with  Target  2  (AF  peak  value 
=  377.45):  (a)  ambiguity  function,  (b)  AF  contour  plot,  (c)  Target  2  frequency 
response  (magnitude),  (d)  transmit  signal  frequency  response  (magnitude),  and  (e) 
return  signal  frequency  response  (magnitude). 


E.  PRACTICAL  TARGET  (BOEING  737-500  RANGE  PROFILE) 

We  now  consider  an  actual  target  response  shown  in  Fig.  5  of  a  Boeing  737-500 
and  we  generate  in  Fig.  11  (for  eigenwaveform)  and  Fig.  12  (for  wideband  waveform) 
the:  a.  AF,  b.  AF  contour,  c.  frequency  response  (magnitude)  of  Boeing  737-500  profde, 
d)  transmit  signal  frequency  response  (magnitude)  and  e.  Boeing  737-500  return  signal 
frequency  response  (magnitude).  The  eigenvalue  of  Boeing  737-500’s  range  profde  in 
Fig.  5  is  60.173  which  leads  to  a  very  high  peak  value  3620.9  in  Fig.  11  while  Fig.  12 
with  wideband  waveform  can  only  achieve  1,  which  is  the  theoretical  result  by  using 
wideband  waveform.  Moreover,  the  eigenwaveform  eliminates  most  Doppler  sidelobes 
but  doubles  the  length  of  return  signal.  In  this  example,  this  type  of  response  is  more  like 
our  Target  2. 
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Boeing  737  Ambiguity  Function 


Boeing  737  frequency  response 


«00- 


•™  -05 


Delays  iti  samples  Doppler  in  Noimal^  Frequency 


Return  signal  frequency  response 

(50 1 - ^ - ; - , - 1 - 1 - ^ 


‘  ‘  1  ' 

_ ^ ^ ^ _ _ 

■05  -or  -03  -02  -or  0  01  02  03  0<  05 


Figure  1 1 .  Ambiguity  function  analysis  of  eigenwaveform  with  Boeing  737-500;  (a) 
ambiguity  function,  (b)  AF  contour  plot,  (c)  Boeing  737-500  frequency  response 
(magnitude),  (d)  transmit  signal  frequency  response  (magnitude),  and  (e)  return 
signal  frequency  response  (magnitude). 
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(a)  Boeing737  Ambiguity  Function  (Wideband) 

1^ 


05-t 


•500  -05 


Delays  in  samples  Doppler  in  Normalized  Frequency 


(c)  Boeing737  frequency  response 


(e)  Return  signal  frequency  response 


Figure  12.  Ambiguity  function  analysis  of  wideband  waveform  with  Boeing  737-500;  (a) 
ambiguity  function,  (b)  AF  contour  plot,  (c)  Boeing  737-500  frequency  response 
(magnitude),  (d)  transmit  signal  frequency  response  (magnitude),  and  (e)  return 
signal  frequency  response  (magnitude). 


In  conclusion  we  notice  that  for  oscillatory  targets  (target  responses  with 
resonances  in  particular  bands)  the  wideband  waveform  tends  to  result  in  a  larger  peak 
than  the  rectangular  waveform.  For  less  oscillatory  targets  in  which  the  particular  Boeing 
737-500  response  is  an  example,  the  rectangular  waveform  tends  to  result  in  a  larger 
peak  than  the  wideband  waveform.  However,  regardless  of  the  target  response,  the 
eigenwaveform  provides  the  largest  AF  peak  compared  to  both  waveforms.  Qualitatively, 
the  overall  AF  sidelobe  suppression  of  the  eigenwaveform  is  clearly  superior  than  both 
waveforms. 
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F. 


AF  ZERO-DELAY  AND  ZERO  DOPPLER  CUTS  FOR  EXTENDED 
TARGETS 


In  ambiguity  function  analysis,  two  cross  sections  are  usually  of  interest:  the  AF 
zero-delay  and  the  zero-Doppler  euts.  The  zero-delay  cut  is  simply  the  AF  eross-seetion 
when  =  0  and  zero-Doppler  eut  is  the  AF  eross  seetion  when  The  cross- 

seetions  are  important  sinee  they  usually  eonvey  delay  (or  range)  resolution  and  Doppler 
(or  veloeity)  resolution.  For  elassieal  waveforms,  the  zero-delay  eut  is  given  by 


(32) 


where 


72^-2* 

N 


and  x[n\  is  the  transmit  signal.  Unfortunately,  this  does  not  work  for  extended  targets. 
However,  reeall  that  we  ineorporated  the  faet  that  the  target  has  finite  extent  in  (16), 
whieh  simplifies  to 


(33) 


where  s{n\  =  x{n\ *  h{n\  is  the  return  echo. 

Notice  that  we  ean  simplify  (33)  to 

IFFT[\s[nf{ , 


dn 


1=1 


(34) 


whieh  states  that  the  AF  zero-delay  is  simply  the  squared-magnitude  of  inverse  (fast) 
Fourier  transform  of  the  squared-magnitude  of  the  transmit  waveform  -  target  response 
eonvolution. 

The  eorresponding  zero-Doppler  cut  is  then  the  squared-autoeorrelation  funetion 
of  the  return  echo  (rather  than  just  the  transmit  signal)  and  is  therefore  given  by 
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|2 


(35) 


n=\ 

The  significance  of  the  AF  zero-delay  and  zero-Doppler  cuts  are  illustrated  using 
transmissions  with  multiple  pulses  which  is  the  topic  of  the  next  section. 
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V.  COHERENT  MULTIPLE-PULSE  TRANSMISSION 


When  multiple-pulse  transmission  or  coherent  pulse  train  is  used,  the  ambiguity 
function  is  the  superposition  of  the  ambiguity  function  of  a  single  transmission  with  some 
amplitude  re-scaling.  From  [12],  the  ambiguity  function  for  a  pulse  train  return  is  given 
by 


1 

L 


£  \xx{^-qTRJd)\ 


q=-(L-\) 


sin[;r/^(L-|g|)7;] 


(36) 


I  |2 

where  \Xi\  is  the  ambiguity  function  of  a  single  transmit  waveform-target  response  pair. 
L  is  the  number  of  pulses  in  a  transmission  and  7]^  is  the  separation  between  pulses  or 
pulse  repetition  interval  (PRI).  The  pulse  repetition  frequency  (PRF)  is  defined  as 

/,=T.  (37) 


In  pulsed  Doppler  radar  system,  one  of  the  key  design  elements  is  PRF.  In  practice,  the 
definition  of  low,  medium,  and  high  PRF  is  truly  application  driven.  This  is  because  the 
PRI  (or  PRF)  dictates  the  desired  or  specified  unambiguous  range.  The  unambiguous 
range  is  given  by 


cTr  _  c 
2  2PRF  ’ 


(38) 


or  the  unambiguous  delay  is  given  by 


T  1 
2  2PRF  ■ 


(39) 


When  using  coherent  pulses,  the  Doppler  resolution  is  basically  half  the  mainlobe 
of  the  zero-Doppler  cut  of  (36).  It  is  given  by 


I 

It, 


(40) 
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In  practice,  a  simple  rule  of  thumb  is  that  medium  PRF  is  a  deeade  larger  than  low  PRF 
and  the  high  PRF  is  a  deeade  larger  than  medium  PRF  but  sueh  rule  is  easily  broken 
depending  on  application.  For  the  convenient  illustration  of  the  PRF  eoneepts  in  our  work 
when  illuminating  an  extended  target,  we  will  simply  define  the  following:  a.  high  PRF  is 
when  is  four  to  ten  times  T  (target  duration),  b.  medium  PRF  is  when  7)^  is  forty 
times  T ,  and  c.  low  PRF  is  when  is  a  hundred  times  T .  Here  our  definition  of  high 
PRF  is  pretty  high.  For  example,  even  with  F^=10F,  the  unambiguous  range  is 
=Tr  /  2  =  10T  /  2  =  5T .  In  other  words,  the  unambiguous  delay  for  this  PRF  can 

accommodate  about  five  target  responses.  This  ehoice  is  not  motivated  by  maximizing  or 
minimizing  the  unambiguous  range  from  some  speeific  application  but  rather  to  show 
what  happens  when  a  (very)  high  PRF  is  lowered  into  medium  PRF  and  then  lowered 
again  to  a  low  PRF.  In  other  words,  we  vary  the  length  of  Tr  to  investigate  how  the 

different  PRFs  affect  the  ambiguity  function  in  terms  of  how  the  volume  inside  the  AF 
changes  and  how  range  and  Doppler  resolution  may  be  affected  as  the  PRI  or  PRF  is 
changed.  As  pointed  out  earlier,  the  volume  of  the  ambiguity  function  may  be 
manipulated  depending  on  the  waveform  ehoice. 

Another  parameter  of  interest  in  a  pulse  train  is  the  duty  cycle  where 

dt=Y  =  T,^fR-  (41) 

Multiple-pulse  transmissions  eonvolving  with  finite-duration  targets  resulting  in  multiple 
eehoes  are  illustrated  in:  Fig.  13a  using  wideband  waform.  Fig.  13b  using  rectangular 
waveform,  and  Fig.  13o  using  eigenwaveform. 
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Figure  13.  Multiple-pulse  transmission  for  extended  target  setting:  (a)  wideband 
waveform,  (b)  rectangular  waveform,  and  (c)  eigenwaveform. 


In  Fig.  14  and  Fig.  15  using  L  =  3  pulses  and  T^=\24  delay  samples,  AF 

contours  are  shown  using  the  wideband,  high-ToTR  rectangular,  and  eigenwaveform 
pulse  trains  for  Target  1  and  Target  2  respectively.  It  is  clear  all  the  waveforms 
(rectangular  waveform,  wideband  and  eigenwaveform)  are  very  high  PRF  (AF  contours 
are  very  close  in  delay  spread  since  Tj^=4T)  regardless  whether  it  is  Target  1  or  2.  In 
Fig.  14  (with  Target  1),  we  can  deduce  from  the  AF  contour  that  the  very  high-PRF 
rectangular  waveform  yields  a  slightly  cleaner  sidelobe  suppression  response  compared 
to  the  wideband  waveform.  However,  its  AF  peak  (which  is  0.0298)  is  lower  than  the  AF 
wideband  waveform  peak  (which  is  1).  We  can  also  deduce  that  the  wideband  waveform 
has  the  tightest  range  (delay)  response.  The  eigenwaveform  produces  the  largest  peak 
(which  is  33.5),  tightest  Doppler  response,  and  we  can  also  see  that  it  qualitatively 
produces  the  cleanest  overall  sidelobe  suppression  level.  In  Fig.  15  with  Target  2,  all  the 
waveforms  have  low  overall  sidelobe  levels  with  the  eigenwaveform  having  the  tightest 
AF  response  in  Doppler  spread.  As  expected  the  wideband  waveform  has  the  tightest  AF 
response  in  the  range  (delay)  domain  but  it  also  has  the  lowest  peak  (which  is  1).  The 
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rectangular  waveform  has  a  decent  peak  (which  is  232.38)  but  the  eigenwaveform  has  the 
largest  peak  (which  is  377.45). 


(a)Traget  1  contour  plot  (Wicteband  waveform.  L=3.  Tr=4T) 
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(b)Traget  1  contour  plot  (Rectangular  wavtform.  L=3,  Tr=4T.  ToTR=0.S4) 
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(c)  Tragel  1  contour  plot  (Eigenwaveform.  L=3.  Tr=4T) 
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Figure  14.  Target  1  ambiguity  function  contour  plots  comparison  (L  =  3,  =  4T):  (a) 

wideband  waveform,  (b)  rectangular  waveform,  and  (c)  eigenwaveform. 


Delays  in  samples 

Figure  15.  Target  2  ambiguity  function  contour  plots  comparison  (L  =  3,  =  4T):  (a) 

wideband  waveform,  (b)  rectangular  waveform,  and  (c)  eigenwaveform. 


36 


A, 


ZERO-DOPPLER  CUT 


Utilizing  our  high-PRF  definition  7"^  =107,  3 -pulse  transmission,  and  Target  1, 

we  illustrate  the  AF  zero-Doppler  cuts  of  the  wideband  waveform,  rectangular  waveform, 
and  eigenwaveform  in  Fig.  16.  It  is  clear  from  Fig.  16  that  the  range  resolution  of  the 
wideband  waveform  is  slightly  better  than  the  other  two  waveforms  as  was  already 
inferred  from  the  AF  contours  of  Fig.  14.  The  peak  is  clearly  apparent  (and  is  the  largest) 
for  the  eigenwaveform  AF.  These  are  the  advantages  of  looking  at  another  perspective 
via  zero-Doppler  cuts  that  may  not  be  evident  from  the  AF  contours. 


Zero-Doppler  cut  of  Target  1  vs  wideband  waveform  (L=3,  Tr=10T,  High  PRF) 


Dtlsrys  in  s«mpl»s 

Figure  16.  Zero-Doppler  cut  of  Target  1  AF  with  various  waveforms.  Upper  panel; 
wideband  waveform  (L  =  3,  high  PRF),  middle  panel:  rectangular  waveform  (L  = 
3,  ToTR=0.64,  high  PRF),  and  lower  panel;  eigenwaveform  (L  =  3,  high  PRF). 

Moreover,  when  eigenwaveform  is  used,  7)^  =  4r  (in  Fig.  14c  and  Fig.  15c)  is  the 
smallest  value  so  that  the  zero-Doppler  cut  matched  filtered  pulse  returns  do  not  overlap. 
Thus,  Tj^  =  AT  is  the  minimum  PRI  ( )  such  that  the  AF  matched-filtered  return 

echoes  from  the  same  target  do  not  overlap.  If  is  increased  with  the  use  of 
eigenwaveform,  then  the  matched  filtered  pulse  returns  start  to  overlap.  Thus,  the 
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maximum  PRF  allowable  if  no  overlap  in  the  eigenwaveform  AF  zero-Doppler  eut  is 
desired  is  given  by 


f  = - =  — . 

J  R  niSX  rri  A  rr< 


T  4T 

^Rmm 


(42) 


The  PRI  =  4r  translates  to  the  aetual  pulse  returns  being  separated  by  2T 

from  the  same  target.  This  is  not  to  be  eonfused  from  =  \.5T  (mentioned  in  an  earlier 

seetion),  whieh  is  the  minimum  time  separation  required  for  two  targets  sueh  that  pulse 
returns  from  the  two  targets  do  not  overlap  (i.e.,  sueh  that  the  two  targets  ean  be 
resolved). 

It  is  interesting  to  note  in  Fig.  16  that  the  ratio  of  eaeh  first  sidelobe  peak  to 
mainlobe  peak  for  all  three  zero-Doppler  cuts  to  be  0.44  for  L  =  3  scenario.  The  ratio 
0.44  is  the  square  of  (2/3).  The  ratio  of  each  sidelobe  peak  to  mainlobe  peak  can  be 
formulated  for  any  L  as 


R  = 


L-i 

\  L  j 


(43) 


where  is  the  ratio  of  side  peak  amplitude  to  main  peak  amplitude. 


Another  useful  comparison  is  to  illustrate  the  AF  zero-Doppler  cuts  for  various 
PRFs  given  a  particular  waveform.  Utilizing  the  eigenwaveform  and  L  =  3 ,  we  show  the 
zero-Doppler  cuts  for  high  PRF,  medium  PRF,  and  low  PRF  in  Fig.  17.  Just  like  in 
traditional  pulsed-Doppler  waveforms,  the  low  PRF  yields  the  best  unambiguous  range. 
However,  just  like  in  pulsed-Doppler  waveforms,  the  choice  between  low,  medium,  and 
high  PRF  has  an  impact  in  the  Doppler  domain,  which  is  considered  shortly. 
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Zero-Doppler  cut  of  Target  2  vs  eigenwaveform  (L=3,  Tr=10T,  High  PRF) 
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Zero-Doppler  cut  of  Target  2  vs  eigenwaveform  (L=3,  Tr=40T,  Medium  PRF) 
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Zero-Doppler  cut  of  Target  2  vs  eigenwaveform  (L=3,  TfIOOT,  Low  PRF) 
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Figure  17.  Zero-Doppler  cut  of  Target  2  AF  with  normalized  eigenwaveforms  and 
various  PRFs.  Upper  panel:  (L  =  3,  high  PRF),  middle  panel:  (L  =  3,  medium 
PRF),  and  lower  panel:  (L  =  3,  low  PRF). 


Also,  we  explore  the  impact  of  increasing  the  number  of  pulses  ( L  =  1, 3, 5  )  given 
a  fixed  PRF.  Again,  utilizing  the  eigenwaveform,  we  illustrate  the  AF  zero-Doppler  cuts 
for  Target  1  with  high  PRF.  Of  course,  L-\  is  an  one-pulse  waveform  which  is  shown 
to  illustrate  what  is  gained  by  coherent  pulse  integration  when  using  multiple  pulses. 
Since  the  PRF  is  fixed,  the  separation  between  lobes  for  L  =  3  and  L  =  5  are  the  same. 
The  main  gain  of  coherent  integration  is  clearly  the  gain,  which  is  ,  but  at  the  expense 
of  higher  first  sidelobe-to-main  peak  ratio  (0.44  for  L  =  3  and  0.64  for  L  =  5)  as 
predicted  by  (43). 

In  general  (when  no  compression  is  used),  the  range  resolution  depends  on 
resulting  extent  of  the  convolution  of  the  transmit  waveform  and  target  response.  In  our 
work,  the  wideband  waveform  has  a  slightly  better  range  resolution  than  the 
eigenwaveform  or  rectangular  waveform  due  to  its  narrow  time  extent.  For 
eigenwaveform,  there  is  a  lower  limit  of  =  4r  for  matched  filtered  echoes  not  to 
overlap  in  AF’s  range  (delay)  domain.  Increasing  L  increases  coherent  integration  gain 
while  sacrificing  sidelobe  peak  to  main  peak  ratio.  Increasing  7)^  (which  translates  to 
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lowering  of  PRF)  increases  unambiguous  range  and  its  impact  on  Doppler  domain  is 
covered  in  next  section. 


Zero-Doppler  cut  of  Target  1  vs  eigenwaveform  (L=1,  Tr=10T,  High  PRF) 
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Figure  18.  Zero-Doppler  cut  of  Target  1  AF  with  un-normalized  eigenwaveforms  and 
various  L.  Upper  panel;  (L  =  1,  high  PRF),  middle  panel:  (L  =  3,  high  PRF),  and 

lower  panel:  (L  =  5,  high  PRF). 


B.  ZERO-DELAY  CUT 

The  zero-delay  cut  shows  the  AF  Doppler  spread  when  the  delay  is  zero.  At  times, 
the  zero-delay  cut  offers  a  perspective  that  is  not  quite  apparent  from  an  AF  contour. 
From  [5]  and  [12],  we  realize  that  the  width  of  the  main  lobe  is  decreased  (i.e.,  Doppler 
resolution  improvement)  by  increasing  the  number  of  pulses  in  a  transmission  by 
observing  the  AF  zero-delay  cut,  which  is  given  by 


sm[7uf^LT^]  2 
sin[;z-/^rj 


(44) 


Employing  the  eigenwaveform,  we  illustrate  in  Fig.  19  the  AF  zero-delay  cuts  for  Target 
1  and  Target  2  with  T  =  1,3,5  using  medium  PRF  with  =  40r  .  We  include  T  =  1  to 

I  |2 

remind  us  the  AF  /^„]|  )  corresponding  to  T  =  1  serves  as  the  envelope  as  dictated 
by  (33).  For  traditional  waveforms  for  point  targets,  the  zero-delay  cut  is  purely  a 
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function  of  the  transmit  waveform.  From  Fig.  19,  it  is  clear  that  extended  target  AF  zero- 
delay  euts  differ  from  target  to  target.  Notice  that  the  zero-Doppler  eut  is  the  squared 
multiplication  of  L  =  1  zero-Doppler  eut  with  a  sine-train.  The  frequency  separation  of 
the  sine-lobes  is  dictated  by  1  /  and  the  half  width  of  the  main  lobe  is  given  by 

=  (45) 

LIr 

whieh  is  also  eonsidered  as  the  Doppler  resolution.  However,  in  (44)  notice  that  the 
separation  of  the  sine-lobes  is  also  a  funetion  of  Tr  and  L .  In  other  words,  Doppler 
separation  (to  avoid  ambiguity)  is  also  a  funetion  of  Tr  and  L  .  Thus,  by  varying  Tr  and 

L  ,  we  can  improve  (or  degrade)  Doppler  resolution  and  ehange  the  frequency  separation 
of  the  sine-lobes.  For  example,  we  illustrate  in  Fig.  20  the  eigenwaveform  AF  zero-delay 
euts  for  Target  1  and  Target  2  with  T  =  1,3,  and  5  but  now  using  high  PRF  (smaller  Tr  ).  It 

is  elear  that  the  higher  PRF  (whieh  is  lower  Tr)  has  wider  separation  but  poorer 
resolution  since  the  bandwidth  (45)  is  increased  by  deereasing  Tr  .  Thus,  to  maintain  the 
same  bandwidth  or  Doppler  resolution  in  (45)  while  deereasing  Tr  (inereasing  sine-lobe 
separation),  it  is  neeessary  to  increase  L .  In  other  words,  Doppler  resolution  is  dietated 
by  the  length  of  total  transmission  T^  =  LTr  as  observed  earlier. 
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Target  1  zero-delay  cut  comparison  (L=1,3,5,  Tr=40T,  Medium  PRF,  Eigenwaveform) 
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Target  2  zero-delay  cut  comparison  (L=1,3.5,  Tr=40T.  Medium  PRF.  Eigenwaveform) 


Doppier  in  Normalized  Frequency 


Figure 


19.  Zero-delay  cut  comparison  with  single,  three,  and  five  retransmissions  with 

low  PRF  using  the  eigenwaveform. 


Target  1  zero-delay  cut  comparison  (L=1,3,5,  Tr=10T,  High  PRF,  Eigenwaveform) 


Target  2  zero-delay  cut  comparison  (L=1,3,5,  Tr=10T,  High  PRF,  Eigenwaveform) 


Figure  20.  Zero-delay  cut  comparison  with  single,  three,  and  five  retransmissions  with 

high  PRF  using  the  eigenwaveform. 
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C.  RANGE  AND  DOPPLER  RESOLUTION  TRADE-OFF 

It  is  now  clear  there  is  an  inherent  trade-off  between  range  and  Doppler  resolution 
by  varying  the  PRF  (and  L)  of  multiple  pulse  transmission  of  eigenwaveforms  (or  other 
waveforms  for  that  matter)  for  extended  targets.  To  illustrate  the  range  and  Doppler 
trade-off  that  eontinue  to  exist  even  for  multiple-pulse  waveform  for  extended  targets,  we 
employ  a  high  PRF  and  a  low  PRF  eigenwaveform  on  Target  1  and  show  the 
corresponding  zero-Doppler  cuts  in  Fig.  21a  and  Fig.  21c,  respectively,  and  zero-delay 
cuts  in  Fig.  21b  and  Fig.  21d  with  various  L  given  a  constant  transmit  energy  constraint 
and  a  constant  total  length  eonstraint  (i.e.,  and  LT^  is  eonstant).  From  Fig.  21a  and 

Fig.  21c,  it  is  clear  that  the  range  resolution  is  the  same  since  the  pulse  width  remains  the 
same;  the  unambiguous  range  is  worse  for  the  high  PRF  than  low  PRF.  From  Fig.  21b 
and  Fig.  2 Id,  the  Doppler  resolution  is  the  same  since  LT^  is  constant  but  the  sine-lobe 
separation  is  larger  for  the  high  PRF  than  low  PRF. 

In  summary,  range  resolution  and  Doppler  resolution  depend  on  the  type  of 
transmit  waveform  and  transmit  signal  length.  For  multiple  pulse  eigenwaveforms, 
increasing  PRF  (deereasing  )  results  in  increasing  Doppler  sine-lobes  separation  but 
increasing  PRF  results  redueing  the  unambiguous  range.  When  larger  L  is  used,  Doppler 
resolution  may  improve  but  larger  L  results  in  more  sidelobes,  which  results  in 
sacrifieing  sidelobe  peak  to  mainlobe  peak  ratio.  It  is  a  classieal  trade-off  where  a 
speeific  applieation  dictates  the  right  choice  for  a  system.  But  we  should  mention  the  fact 
that  the  narrowband  nature  of  the  eigenwaveform  (as  opposed  to  the  wideband 
waveform)  helps  in  the  reduction  of  the  sine-lobe  peaks. 
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(a)  Zero-Doppler  cut  of  Target  1  (L=30.  Tr=10T,  High  PRF)  (c)  Zero-Doppler  cut  of  Target  1  (L=3.  Tr=100T,  Low  PRF) 


Figure  21 .  PRF  and  pulses  number  eomparison  with  eigenwaveform;  (a)  zero-Doppler 
cut  (L  =  30,  high  PRF),  (b)  zero-delay  cut  (L  =  30,  high  PRF),  (c)  zero-Doppler 
cut  (L  =  3,  low  PRF),  and  (d)  zero-delay  cut  (L  =  3,  low  PRF). 
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VI.  DETECTION  PROBABILITY  AND  RANGE-DOPPLER  MAP 


In  this  section  we  investigate  the  probability  of  detection  of  various  waveforms 
such  eigenwaveform,  rectangular  waveform  and  wideband  waveform.  We  can  assume  the 
target  to  have  zero  Doppler  and  its  extension  to  the  generalized  case  is  straightforward. 


A,  DETECTION  PROBABILITY  OF  BASIC  ONE-PULSE  WAVEFORMS 


Let  h  be  the  complex -valued  target  response  and  w  be  the  complex  valued  white 
Gaussian  noise  in  the  receiver  with  a  sample  variance  of  cr^ .  Let  s  be  the  convolution  of 
transmit  waveform  x  and  h .  Then  the  detection  hypotheses  are  given  by 


^o:y  =  w 

//j :  y  =  s  +  w  =  Hx  +  w 


(46) 


where  H  is  the  target  convolution  matrix  corresponding  to  h.  The  decision  statistic 
using  matched  filter  theory  for  a  fixed  threshold  y  is 

r(y)  =  Re{y"s}  =  Re{y"Hx}.  (47) 

When  the  wideband  waveform  is  used,  the  detection  probability  [18]  given  a  fixed  false 
alarm  probability  (PFA)  is  given  by 

J2F  F 

(48) 

=  Q{Q-'(P,J-pE,mR)),  (49) 


where  TNR  =  /  cr^  (received  target-to-noise  ratio)  and  Q  function  is  the  tail 

probability  of  the  standard  normal  distribution.  Alternatively,  (22)  can  be  given  as 

P,=Q{Q-'{P„)-pE,SNR)),  (50) 


where  SNR  =  (transmit  signal-to-noise  ratio).  In  other  words,  is  a  function  of 

both  transmit  energy  and  received  target  energy  along  with  the  noise  variance.  A  more 
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compact  version  of  (48),  which  resembles  detection  probability  of  traditional  matched 
filter  analysis  [19]  is 


Po=Q(Q^'(P,J-pB,NR)l 


(51) 


where  E^NR  =  /  cr^  where  specifically  means  return  (eeho)  energy  using  the 

wideband  waveform.  It  can  be  shown  [7]  that  the  detection  probability  with  the  use  of 
eigenwaveform  is  given  by 

P^=Pr{T>r;H,) 

= = Q(  ^  ) 

V  2  V  2  (52) 

Notice  that  performance  improvement  with  the  use  of  eigenwaveform  compared  to  the 
wideband  (impulse)  waveform  given  a  fixed  PEA  by  comparing  (25)  and  (26)  where 

effectively  amplifies  the  E^NR  in  (51). 


For  any  arbitrary  waveform  (rectangular  waveform  ineluded),  it  can  be  shown  that 
the  generalized  detection  probability  is  given  by 


F„=Pr(r>y;//,) 


icT^. 

y 


—H  —  — 

2^Ay 

X  RaX 

(53) 


where  x  and  Ra  are  normalized  transmit  waveform  and  target  response  autocorrelation 


matrix.  It  is  intuitive  to  conelude  that  the  value  of  ( 


—H —  — 

X  Rax 


)  may  vary  from  close  to  zero 


(since  Ra  is  positive  definite)  to  As  a  consequence,  the  detection  performance 
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pretty  much  depends  on  the  how  good/bad  the  transmit  waveform  matches  target 
response’s  autocorrelation  matrix.  Thus,  if  we  want  to  compare  detection  probabilities  of 
arbitrary  waveforms,  they  have  to  be  calculated  via  (53)  and  only  then  we  can  tell  if  it  is 
greater  or  less  than  the  detection  probability  of  the  wideband  impulse  waveform  that  is 
given  in  (48).  In  other  words,  stating  what  the  transmit  energy  is  not  enough  to  know  the 
detection  probability  of  an  arbitrary  waveform.  Its  correlation  with  the  target  dictates  the 
total  return  energy.  The  eigenwaveform  ensures  the  maximum  detection  probability  since 

—H —  — 

the  term  x  Rax  inside  the  square  root  in  (53)  is  maximized  which  yields  the  detection 
probability  in  (52).  The  detection  performance  curves  for  various  waveforms  illuminating 
Target  1  and  Target  2  are  shown  in  Fig.  22  and  Fig.  23  given  a  fixed  TNR  (target-to- 
noise  energy  ratio)  for  various  false  alarm  probabilities  Considering  Fig.  22,  the 

performance  of  the  eigenwaveform  is  superior  to  rectangular  and  wideband  waveforms 
(given  a  fixed  as  expected.  For  example,  P^^=0.01  and  P^=Q.9,  the  eigenwaveform 
advantage  over  the  wideband  waveform  is  about  7.6  dB  (which  makes  sense  since  this  is 
about  equal  to  101og(/f^jj^).  Whaf  s  more  incredible  however  (for  oscillatory  Target  1)  is 

that  the  eigenwaveform  with  a  stringent  requirement  of  P^^  =  1 0  still  outperforms  both 

the  wideband  waveform  and  the  rectangular  waveform  with  a  looser  requirement  of  1 0^^ 
(for  P^>0.1)!  In  Fig.  23  where  Target  2  is  considered,  the  performance  curve  for 
eigenwaveform  still  is  the  best  compared  to  rectangular  waveform  and  wideband 
waveform  given  a  fixed  P^^ .  However,  as  can  be  inferred  from  the  ambiguity  functions 

previously  studied,  the  detection  performance  of  the  rectangular  waveform  is  much  better 
than  the  wideband  waveform  (for  the  non-oscillatory  Target  2).  This  reinforces  the  notion 
that  was  stated  earlier.  We  know  that  the  eigenwaveform  yields  the  best  detection 
performance  but  detection  performance  of  other  waveforms  have  to  be  calculated  via  (53) 
such  that  performance  comparison  (in  terms  of  detection)  between  waveforms  can  be 
made.  Finally,  the  detection  performance  comparison  adds  the  needed  final  dimension  to 
the  range -Doppler  trade-off  when  considering  various  waveforms.  In  other  words, 
systems  which  are  noise-limited  may  opt  for  the  eigenwaveform  with  a  slight  hit  on  range 
resolution.  Other  systems  which  may  not  be  noise-limited  may  opt  for  other  waveforms  if 
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range  resolution  takes  precedence  (but  with  possibly  substantial  cost  in  detection 
performance). 
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Figure  22.  Target  1  detection  probability  comparison  of  wideband  waveform,  rectangular 
waveform  (ToTR  =  0.64)  and  eigenwaveform  (TNR  =  OdB). 


Target  2  (Non-oscillatory  target)  Detection  Performance  Comparison; 
Eigenwaveform,  Rectangular  (ToTR=0.64)  and  Wideband  waveform  (TNR  OdB) 
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Figure  23.  Target  2  detection  probability  comparison  of  wideband  waveform,  rectangular 
waveform  (ToTR  =  0.64)  and  eigenwaveform  (TNR  =  OdB). 
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B,  DETECTION  PROBABILITY  OF  COHERENT  MULTIPLE-PULSE 

TRANSMISSION  AND  RANGE  DOPPLER  MAP  APPLICATION 

When  moving  targets  are  present,  one  of  the  ambiguity  function  application  is  the 
range-Doppler  map  (RDM)  which  significantly  demonstrates  the  benefit  of  using 
eigenwaveform  for  extended  targets.  Assume  L  pulses  are  sent  and  that  the  return  echo 
is  received.  Any  target  detected  may  be  moving  such  that  a  Doppler  component  is 
possible.  After  receiving  the  return,  the  long  sequence  is  carefully  re-arranged  so  that 
each  return  from  every  pulse  is  aligned  according  to  same  delay  and  stored  in  a 
measurement  matrix.  By  taking  fast-Fourier  transform  (FFT)  in  Doppler  direction,  signal 
energy  converges  into  corresponding  Doppler  bins  (indices)  due  to  FFT’s  circular  shift 
property.  The  magnitude  of  the  measurement  matrix  after  the  FFT  operation  is  considered 
the  range  Doppler  map  since  it  reveals  the  characteristic  of  moving  target’s  delay  (which 
corresponds  to  range)  and  Doppler  shift  (which  corresponds  to  velocity).  The  multiple- 
pulse  wideband  waveform,  rectangular  waveform,  and  eigenwaveform  are  utilized  as  the 
transmit  waveforms  here.  A  peak  value  may  be  used  to  detect  a  target’s  position  (delay) 
and  Doppler  spread.  Of  course,  the  RDM  can  detect  multiple  targets  with  different  speeds 
and  distances  and  ambiguities  can  be  avoided  via  discussion  from  earlier  sections. 

Moreover,  when  multiple  pulses  are  used,  it  can  be  shown  that  the  detection 
probability  is  also  a  function  of  the  number  pulses  (L  ).  The  detection  probability  with  L 
pulses  with  the  use  of  eigenwaveform  is  be  shown  to  be 

P^=Vx{T>r-H,) 

= 

V  2  V  2 

The  gain  due  to  L  (10  logL  in  dB)  is  also  true  for  other  waveforms  as  shown  in  Fig.  24. 
In  Fig.  24,  we  illustrate  the  Target  1  detection  performance  curves  for  the  three 
waveforms  as  a  function  of  increasing  L  .  As  expected,  increasing  L  increases  detection 
probability. 


)  =  Q{Q\p,a) 


21 


cr^ 


(54) 
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Now,  we  generate  the  RDM  (3-D  maps)  where  we  utilize  Target  1  and  illuminate 
it  with  the  wideband,  reetangular  (ToTR  =  0.64),  and  eigenwaveform  (duty  eycle  d,  = 

10^^)  where  Z  =  31  as  a  funetion  of  inereasing  SNR  (0,  3,  10  dB)  given  a  fixed  TNR  =  0 
dB.  We  illustrate  in  Fig.  25  the  3-D  RDMs  of  these  various  waveforms  as  a  function  of 
increasing  (return  energy  to  noise  ratio).  In  this  scenario,  the  target  is  located  at 

range  (delay)  =  560  with  normalized  Doppler  of  =  0.0968 .  Notice  that  the  target  is  not 
very  visible  for  SNR  =  0  dB  with  the  use  of  rectangular  waveform.  For  the  same  SNR, 
there  is  a  small  peak  corresponding  to  Target  1  with  the  use  of  the  wideband  waveform. 
Notice  however  that  the  Target  1  is  clearly  discernable  with  SNR=0  dB  when  the 
eigenwaveform  is  used.  Increasing  SNR  enhances  the  peaks  for  all  waveforms  with  the 
eigenwaveform  clearly  yielding  the  largest  peak. 

RDMs  are  usually  presented  as  2-D  maps  (magnitude  or  squared-magnitude).  In 
this  scenario,  we  place  two  targets  (Target  1  type)  in  two  different  range-Doppler  bins. 
Again  we  use  the  wideband,  rectangular  (ToTR  =  0.64),  and  eigenwaveform  (t/,  =10^^) 
where  Z  =  31  as  a  function  of  increasing  SNR  (0,  3,  10  dB)  given  a  fixed  TNR  =  OdB. 
The  RDMs  are  illustrated  in  Fig.  26.  In  this  scenario,  the  targets  are  located  at  range 
(delay)  =  450  with  normalized  Doppler  of  (/^j  =0.1613)  and  range  (delay)  =  750  with 
normalized  Doppler  of  (/^j  =  “0-1290).  Again,  notice  that  the  targets  are  not  very 

discernable  for  SNR  =  0  dB  with  the  use  of  rectangular  waveform.  For  the  same  SNR, 
there  is  a  some  faint  possibility  that  one  of  the  targets  maybe  detected  with  the  use 
wideband  waveform.  Both  targets  are  clearly  pronounced  when  the  eigenwaveform  is 
used  even  with  SNR  =  0  dB.  Increasing  SNR  enhances  the  target  bins  in  the  RDMs  for  all 
waveforms  with  the  eigenwaveform  clearly  yielding  the  brightest  target  bins. 
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Target  1  Detection  Performance  Comparison:  (TNR  OdB,  PFA=10'®) 


Figure  24.  Performance  comparison:  wideband  waveform,  rectangular  waveform  and 
eigenwaveform  with  various  multiple-pulses  transmission  (L  =  3,  5,  10). 
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Figure  25.  Target  1  3-D  view  RDM  comparison  for  SNR  =  (0  dB,  3  dB,  10  dB),  and 
TNR  =  0  dB.  The  target  is  located  at  range  (delay)  =  560  with  Doppler  shift  of 

=  0.0968.  Left  panels:  wideband  waveform,  middle  panels:  rectangular  waveform, 

and  right  panels:  eigenwaveform. 
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Figure  26.  Target  1  RDM  comparison  for  SNR  =  (0  dB,  3  dB,  10  dB),  and  TNR  =  0  dB. 
The  targets  are  located  at  range  (delay)  =  450  with  Doppler  shift  of  /^j  =  0.1613 

and  range  (delay)  =  750  with  Doppler  shift  of  =  0.1290.  Left  panels; 

wideband  waveform,  middle  panels;  rectangular  waveform,  and  right  panels; 

eigenwaveform. 
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VII.  PROBABILITY  WEIGHTED  EIGENWAVEFORMS  (PWE) 
FOR  TARGET  RECOGNITION 


A  closed-loop  radar  which  uses  adaptive  waveforms  is  used  in  [2,  8]  for  target 
recognition.  Among  these  waveforms,  the  PWE  scheme  introduced  in  [1]  has  shown  the 
best  promise  as  far  as  classification  performance  as  well  as  low  mean  number  of 
transmissions.  Here  we  investigate  maximum  a  posted  PWE  (MAP-PWE).  Based  on  the 
decision  statistics,  we  propose  an  improved  un-biased  version  called  match-filtered  PWE 
(ME-PWE).  Moreover,  we  also  propose  two-stage  PWE  (TS-PWE)  which  is  based  on 
likelihood  and  waveform  weighting  that  effectively  works  well  in  low  SNR.  We  describe 
the  closed-loop  radar  operation  involving  these  waveforms  and  summarize  the  radar 
operational  steps  in  a  flow  diagram.  We  perform  Monte  Carlo  simulations  to  show  the 
identification  performance  comparisons. 

A.  ADAPTIVE  WAVEFORMS,  LIKELIHOOD  UPDATES,  AND  COGNITIVE 

RADAR  PLATFORM  INTRODUCTION 

What  makes  a  radar  closed-loop  is  the  use  of  probability  or  likelihood  updates 
into  forming  adaptive  waveforms.  Here  we  consider  the  use  of  closed  loop  or  cognitive 
radar  platform  for  target  recognition  which  had  been  extensively  used  in  [1,  2,  and  8]. 
Since  the  notion  of  adaptive  waveforms,  probability/likelihood  updates  and  waveform 
weighting  are  so  integrated  in  the  closed-loop  radar,  we  will  introduce  our  first  adaptive 
waveform  MAP-PWE  and  the  cognitive  radar  operations  simultaneously  in  the  following 
section. 


1,  MAP-PWE 

Consider  a  target  identification  problem  in  which  one  of  M  possible  targets  is 
present.  Each  target  hypothesis  is  characterized  by  its  impulse  response  = 
which  are  assumed  a  priori. 

Assuming  all  targets  are  of  length  n  and  assuming  there  is  one  target  present,  the 
detection  hypotheses  are 
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//j:y  =  x*hj+w  =  HjX  +  w 
//,  :y  =  x*h,  +  w  =  H,x  +  w 

"  '  (55) 


Hm  :y  =  x*h^+w  =  H^x  +  w 


where  the  H^.  is  the  2n-by-n  convolution  matrix  of  target  j  and  w  is  the  complex¬ 
valued  AWGN  vector  with  a  sample  variance  of  cr^ .  The  correlation  matrix  of  w  is 
simply  cr^I  and  we  assume  cr  to  be  one  for  simplicity.  Each  time  the  radar  transmits  a 
waveform  x ,  a  noise-corrupted  version  of  the  reflected  target  echo  is  received  as  y . 

In  PWE,  the  transmit  signal  is  the  linear  combination  of  each  unit  energy 
eigenvector  q .  each  weighted  by  where  w  .  is  the  weight  distribution  for  the  j“' 

hypothesis  calculated  from  prior  received  signal.  The  energy  constraint  of  transmit  signal 
in  the  end  actually  dictates  the  energy  weight  distribution  of  each  eigenwaveform. 

However,  to  simplify  the  procedure  it  is  best  to  simply  add  the  weights  to  sum  up  to  1  (as 
an  intermediate  normalization  step).  In  other  words, 

M 

2;>v,=i,  (56) 

7=1 

and 


X  = 


(57) 


where  x  is  the  direct  combination  result  of  all  eigenvector  q^  .  However,  x  may  not  be 

unit-energy  after  summation  since  the  eigenwaveforms  themselves  come  from  different 
targets  and  may  not  be  orthogonal.  Although  the  energy  may  be  close  to  1,  it  still  needs  to 
be  normalized  such  that  the  energy  constraint  can  easily  be  factored  in  i.e.. 


X  = 


(58) 


where  is  the  energy  of  transmit  signal  and  E-  is  the  energy  of  x  in  (57). 
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Now,  we  discuss  how  the  waveform  weights  are  calculated.  Let  /;|,  (y)  be  the 

probability  density  or  likelihood  function  of  the  f'  hypothesis  given  i‘^  target  is  the 
present  target  assuming  additive  white  Gaussian  noise.  Recalling  <7  =  1  for  convenience, 
the  probability  density  or  likelihood  function  is  given  by 

fj\i  (y)  =  P  exp[-(y  -  )"  (y  -  s^. )] 

=  /?exp[-y"y  +  s"y  +  y"s^.-s"sJ  (59) 

=  P  exp[-y"y  +  2  Re(s"y)  -  s"s  J, 

where 

(M) 

n 

is  the  constant  in  front  of  the  Gaussian  distribution. 

Let  w]  be  the  initial  waveform  weight  for  each  hypothesis  as  dictated  by  (57).  If 
there  is  no  a  priori  information  available  as  to  the  likelihood  of  each  hypothesis,  then 
initially  we  can  assume  them  to  be  equally  likely  i.e.,  w\  =1!  M .  Let  (y)  be  the 

likelihood  function  from  P‘''  return  signal,  then  the  waveform  weights  are  updated 

by  the  likelihood  values  from  the  latest  R'*  return  signal.  In  other  words,  for  multiple 
illumination  the  waveform  weights  are  updated  as: 

1  1 

w.  =  — 

'  M 

=  cc'wj! 

(61) 

P+\  P  P  rP 

where  fj’  is  the  calculated  likelihood  value  after  P‘^  transmission  and  the  weight  is 
the  weight  distribution  corresponding  to  the  j"‘  hypothesis  for  the  (R  +  1)'*  transmission 
(or  P  updates)  while  ensures  unity  weight  summation  as  dictated  by  (56)  in  each 
transmission.  In  other  words,  is  the  waveform  weight  of  f*'  eigenwaveform  in  (57). 
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It  can  be  shown  that  by  substituting  y  =  H.x  +  w  and  s^.  =H^.x  into  (59),  the 
generalized  formula  of  the  likelihood  function  for  f’'  hypothesis  given  target  is 

(x)  =  P  exp[-x"Hf  H .X  -  x"Hf  w  -  w^H .x  -  w^w  +  2  Re(x"H"H .x  +  x"H" w) 


-x"H"H  X], 

and  the  first  moment  (mean)  of  incorrect  hypotheses  (where  j)  is 


(x)]  =  exp[;7  +  2  Re(x^H;H,x)  -  x" R  .x] 


(62) 


(63) 


where  is  the  discrete  autoeorrelation  function  corresponding  to  j'  hypothesis  and  rj 
is  defined  as 


T]  =  -x“H“H.x-£'[w"w] 


(64) 


whieh  is  a  constant  term  for  all  hypotheses. 

On  the  other  hand,  the  first  moment  (mean)  of  correct  hypothesis  where  i  =  j 
becomes 


EUi,  (X)]  =  /?  exp[/7  +  2  Re(x"Hf  H,x)  -  x^Hf  H,x] 
=  /?exp[77  +  2Re(x"R.x)-x"R.x]. 


(65) 


It  is  clear  from  (65)  that  the  argument  inside  the  exponential  of  the  correct  hypothesis 
{i  =  j)  is  greater  than  all  incorrect  hypotheses  in  (63)  sinee  the  second  term  in  (65)  is  a 

fully-correlated  match  (autocorrelation  matrix  HfH;)  while  (63)  involves  cross-target 
correlation  which  is  clearly  drives  it  to  a  lower  value.  Thus,  the  likelihood  of  eorrect 
hypothesis  has  higher  value  whieh  leads  to  higher  weight  of  w.  in  next  waveform 

formation  in  (61).  Here  we  use  w,.  since  j  =  i  for  the  correct  hypothesis.  Therefore,  the 
waveform  weight  updating  procedure  leads  to  likely  identification  until  a  desired 
probability  threshold  is  met  or  at  the  end  of  desired  number  of  transmissions.  This  is  how 
the  MAP-PWE  cognitively  approaches  the  eorrect  hypothesis  during  its  operation. 
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MAP-PWE  has  been  shown  to  be  a  very  effeetive  adaptive  waveform  formation 
seheme  for  non-moving  target  identifieation  in  this  elosed-loop  radar  platform.  However, 
when  eomparing  (63)  and  (65),  the  last  terms  cause  undesirable  effect  of  subtracting  a 
higher  value  in  the  correct  hypothesis  than  in  the  incorrect  hypotheses.  To  illustrate  this 
effect,  consider  the  case  when  the  updating  procedure  finally  homes  in  on  the  target 
present  (after  a  few  iterations).  In  other  words,  when  w.  « 1  and  x  «  ,  the  last  term 

of  incorrect  y'*  hypothesis  in  (63)  is  an  insignificant  value  since  x  =  is  not  highly 
correlated  to  any  target  j  while  correct  hypothesis  subtracts  a  maximum  value  of 
/l-max  i^x  •  The  biased  subtraction  of  this  term  in  the  most  likely  hypothesis  affects  the 
identification  performance. 

To  illustrate  the  biased  effect  in  MAP-PWE,  consider  the  noise  only  scenario 
where  y  =  w  and  =  H^x  .  It  can  be  shown  from  (59)  that  the  first  moment  of  likelihood 

function  for  all  hypotheses  fj  (x)  is 

E{f  .  (y)]  =  y5exp[-E'[w^w]-x"H"H^x]  =  y5exp[-E'[w^w]-x"R^.x].  (66) 


Since  x  is  the  combination  of  all  eigenwaveforms  with  different  weight  ,  we  can 
further  derive  (66)  in  terms  of  eigenvalues  and  waveform  weights,  that  is 


E[fj  (y)]  =  exp[-E[w'"w]  -  ^ 

Ex 


(67) 


Due  to  the  fully-correlated  match  (autocorrelation  R^. )  inside  the  argument  of  the 
exponential  and  after  a  few  iterations,  the  waveform  weight  updating  procedure 
approaches  the  target  hypothesis  of  higher  eigenvalue  X. .  The  y'*  maximum 

eigenvector  ({  ■  dominates  the  f  j  (x) ,  that  is 


EUi  (y)]  =  y0exp[-E[wV]-|^w^.qfR,qJ 

=  p exp[-E[w'^w] - ^ X.(\  . ]  =  p exp[-E[w"w] - ^ w.Xjl 


(68) 
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Since  the  mean  in  (68)  is  different  for  each  hypothesis,  the  waveform  weights  dictated  by 
(61)  and  (68)  will  have  unequal  weights  after  each  update.  For  the  noise-only  ease 
scenario  that  we  are  currently  discussing,  no  eigenwaveform  should  be  favored  since  a 
target  is  not  present.  Thus,  the  probability  weight  updating  procedure  under  noise  only 
scenario  is  biased  (a  function  of  target  response  matrix’s  eigenvalue)  and  the 
identification  performance  when  a  target  is  truly  present  suffers  from  the  various 
eigenvalues  in  (68).  In  fact,  the  biasing  effect  becomes  worse  when  energy  of  x  increases 
(i.e.,  when  the  transmit  energy  of  the  PWE  waveform  is  increased).  To  illustrate  this 
point,  we  perform  Monte  Carlo  (MC)  trials  where  there  are  four  target  hypotheses  where 
we  assume  initial  probability  of  0.25  for  each  hypothesis  but  under  the  noise  only 
scenario.  Since  there  is  not  a  target  present,  we  only  receive  noise  in  the  measured  signal 
and  calculate  the  4  likelihood  and  waveform  weights.  We  perform  10000  MC  trials  in 
which  to  average  over.  In  Fig.  27a,  we  show  the  waveform  weight  updates  as  a  function 
of  increasing  transmission  energy  using  MAP -PWE.  At  low  transmit  energy,  the  average 
updated  waveform  weight  remains  at  0.25  which  is  obviously  desired.  Notice  however 
that  as  transmission  energy  is  increased,  the  updated  weights  diverge  where  one 
hypothesis  seems  to  be  favored  than  others  (which  is  clearly  unwanted). 

(a)  Mean  of  probability  updates  with  MAP-PWE  under  the  noise-only  scenario 


(b)  Mean  of  probability  updates  with  MF-PWE  under  the  noise-only  scenario 


Energy  per  transmission  (dB) 


Eigure  27.  Mean  of  target  probability  updates  versus  transmission  energy  under  noise 
only  scenario  assuming  unity  noise  energy.  (a)MAP-PWE.  (b)ME-PWE. 
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2. 


MF-PWE 


It  is  our  desire  to  improve  MAP-PWE.  Clearly,  one  way  is  to  remove  the 
unnecessary  bias  terms  in  the  likelihood  functions.  In  reality,  removing  bias  terms  lead  to 
a  receiver  configuration  that  is  equivalent  to  using  a  bank  of  matched  filters.  Removing 
the  constant  terms  in  (59)  lead  to  a  modified  likelihood  function  given  by 

fj\i  (y)  =  P  exp[2  Re(s"y)],  (69) 

where  (5  now  here  is  a  constant  that  ensures  the  waveform  weights  sum  up  to  1  as 
dictated  by  (56).  By  substituting  y  =  H.x  +  w  and  s^.  =H^.x  into  (69)  where  it  is  clear 
why  the  current  technique  is  to  be  called  matched-filter  PWE  (ME-PWE),  (69)  becomes 

//i,'  (y)  =  y0exp[2Re(x"HfH.x  +  x'^Hfw)]. 

It  can  be  shown  the  mean  of  likelihood  for  the  incorrect  hypotheses  given  target  is 
given  by 

IF  {  M  M  ) 

n/ji,  (y)]  =  /9exp[2Re(s''y)]  =  /9exp[^Re  (£»>,q  ,)“H“H,(£».,q,)  ],  (70) 

y  y=i  y=i  J 

and  the  mean  of  the  likelihood  for  the  correct  hypothesis  probability 

2.E  {  M  M  1 

n4  (y)]  =  /»exp[2Re(s''y)]  =  /Sexp[^Re  (y;w,q,)''R,(y;w,q,)  ].  (71) 

y  j=\  7=1  j 

Notice  that  from  (71)  that  the  argument  of  the  exponential  of  correct  hypothesis  {i  =  y  )  is 
ensured  to  be  greater  than  all  incorrect  hypotheses  in  (70)  since  the  fully-correlated  match 
term  in  (71)  surpasses  the  cross-target  correlation  matrix  H^H;  in  (70).  Thus,  the  higher 

value  of  likelihood  for  the  correct  hypothesis  leads  to  higher  waveform  weight  towards 
the  correct  target. 

Moreover,  after  few  updates,  when  the  weight  is  mainly  distributed  to  the 
eigenwaveform  corresponding  to  the  correct  hypothesis,  the  expression  in  (71)  is 
dominated  by 
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(72) 


EUi,  (y)]  -  y5exp[^(qf  R,q,.)]  =  ^exp[^M], 

h-  h- 

X  X 

and  the  means  of  the  likelihoods  of  the  ineorrect  hypotheses  are  given  by 

IF  c  ^ 

£[/,i,(y)]  =  /Sexp[2Re(s"y)]«/5exp[-^Re  q,''HfH,q,  ], 

whieh  leads  to  a  mueh  smaller  value  than  (72). 

We  now  eompare  MAP-PWE  and  MF-PWE  under  noise-only  seenario.  The 
modified  likelihood  f.  (y)  is 

/.(y)  =  ;0exp[2Re(s;y)]  =  y0exp[2Re(x"H»]  (73) 

and  the  mean  is 

EUi  (y)]  =  P  exp[2  Re(x"Hf  E[w])]  =  fi.  (74) 

Clearly,  the  means  of  the  likelihood  functions  are  the  same  under  the  noise-only 
scenario.  Since  we  sum  up  the  likelihoods  and  normalize  to  1  to  produce  the  waveform 
weights,  each  eigenwaveform  will  have  equal  weights  i.e.,  no  target  is  favored  when  a 
target  is  not  present. 

Eastly,  the  waveform  weight  update  rule  remains  the  same,  i.e., 

1  1 
w.  = - 

'  M 

2  I  \  r\ 

W.  =awj^ 

3  2  2  c\ 

Wi  =  a  wj^  (75) 

P+l  P  P  rP 

w,  =  «  w,  / 

where  fP  is  the  calculated  likelihood  weight  corresponding  to  the  P  hypothesis  for  the 

P"'  transmission  while  ensures  unity  weight  summation  in  each  transmission. 

Comparing  (68)  with  (74)  for  noise-only  scenario,  it  is  obvious  that  MF-PWE 
provides  constant  means  for  all  hypotheses  no  matter  how  large  the  eigenvalues  are. 
Thus,  the  MF-PWE  is  an  unbiased  target  identification  scheme  for  extended  targets 
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where  only  the  correet  target  ean  make  significant  likelihood  increase  without  the 
degrading  bias  terms.  Also,  MF-PWE  reduces  the  complexity  of  calculation  used  by 
MAP-PWE  by  simply  using  the  matched  filter  bank  in  (69). 

The  performance  comparison  of  MAP-PWE  and  ME-PWE  vs  various 
transmission  numbers  are  shown  in  Fig.  28  and  Fig.  29  via  10^  Monte  Carlo  experiments 
where  one  true  target  is  present  from  four  hypotheses.  When  the  number  of  transmissions 
is  fixed,  the  hypothesis  with  the  largest  updated  likelihood  is  decided  to  be  the  correct 
hypothesis  (whether  true  or  not).  It  is  clear  in  Fig.  28  that  MF-PWE  has  the  better 
identification  performance  than  MAP-PWE  for  various  numbers  of  fixed  transmissions.  It 
is  also  clear  in  Fig.  29  that  MF-PWE  has  lower  required  mean  number  of  transmissions  to 
reach  a  desired  probability  threshold  (for  Fig.  29  the  threshold  is  set  as  0.95).  In 
conclusion,  MF-PWE  performs  better  than  MAP-PWE  and  wideband  waveform  as  a  CR 
adaptive  waveform  in  terms  of  target  identification  and  required  number  of  transmissions. 


MAP-PWE  vs  MF-PWE  performance  with  multiple  transmissions 


Eigure  28.  Performance  comparison;  MAP-PWE  and  MF-PWE  with  multiple 

transmissions  assuming  unity  noise  energy. 
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PWE  Average  number  of  iterations  vs  transmit  energy  performance 


Figure  29.  Average  or  mean  number  of  transmit  (probability  threshold  =  0.95) 

comparison  assuming  unity  noise  energy. 

B.  TWO-STAGE  PWE  FOR  LOW  SNR 

The  MAP -PWE  and  MF-PWE  closed-loop  radar  adaptive  waveform  techniques 
operates  on  the  previous  likelihood  and  latest  likelihood  update  which  is  primarily  based 
on  Bayes’  theorem.  We  propose  a  new  update  process  that  uses  the  ratio  of  the  latest 
calculated  likelihood  to  the  previous-to-last  calculated  likelihood.  In  low  SNR,  the  noise 
realizations  are  large  (compared  to  actual  target  return)  which  may  cause  an  incorrect 
hypothesis  to  have  a  large  likelihood  update.  When  an  incorrect  hypothesis  takes  on  a 
large  likelihood  update,  it  is  hard  to  reduce  the  next  waveform  weight  since  the  previous 
update  is  incorporated  into  the  actual  waveform  weight.  In  other  words,  if  a  receiver  were 
to  stop  to  make  a  decision  and  the  incorrect  hypothesis  happens  to  have  the  largest 
likelihood,  then  the  receiver  will  choose  that  incorrect  hypothesis.  The  key  idea  is  to 
mitigate  large  swings  on  the  likelihood  updates  (in  the  case  of  low  SNRs)  where  large 
noise  realizations  produce  volatility  to  the  waveform  weights  (which  can  nullify  the 
beneficial  improvements  of  multiple  transmissions). 
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In  this  section,  we  propose  a  heuristic  two-stage  PWE  method  (TS-PWE)  which 
is  specifically  designed  to  improve  the  classification  performance  over  the  MAP-PWE 
and  ME -PWE  in  low  SNR  by  using  the  previous-to-last  and  latest  likelihood  ratios  in  a 
specific  way.  To  describe  the  waveform  technique,  let  us  briefly  go  back  the  closed-loop 
radar  system  that  we  have  been  using  and  see  how  the  two-stage  PWE  is  incorporated. 
The  steps  are: 

1.  Assuming  M  possible  extended  targets,  generate  the  normalized 
eigenwaveform  for  each  target  hypothesis.  Scale  each  eigenwaveform  with  the  square- 
root  of  the  initial  waveform  weight  assigned  to  each  hypothesis  and  then  form  the  first 
PWE -based  waveform. 

2.  After  illumination  and  reception  of  target  return,  calculate  the  likelihood 
according  to  (69).  If  this  is  the  first  transmission,  then  calculate  the  second  set  of 
waveform  weights  via  the  wf  in  (75).  Our  goal  is  to  use  the  latest  and  previous  weight 
pairs  to  generate  the  next  set  of  waveform  weights  for  the  next  transmission. 

3.  Eor  the  second  and  all  the  rest  transmissions,  upon  the  reception  of  target 
return,  pick  the  hypothesis  with  the  highest  weight  (from  the  latest  or  last  weight  set)  and 
calculate  the  ratio  of  weight  change  of  this  very  hypothesis  as  given  by 

where  p  denotes  the  latest  weight  and  p-\  denotes  the  weight  previous  to  last. 

4.  Pick  the  latest  and  p-\  (previous-to-last)  likelihoods  from  the  same 
hypothesis  and  calculate  the  likelihood  ratio  as  given  by: 

Hp) 

y  —  _ 

where  /  is  from  (69). 

5.  Multiply  the  latest  by  the  corrective  coefficient  which  is  the  ratio  of 
over  r,  to  form  the  tempered  or  modified  likelihood  as  given  by 
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y(p)  _  y(p)  Zk 

n 

6.  Calculate  the  new  weight  using  (61)  with  the  modified  likelihood.  Since  the 
weights  sum  up  to  1,  the  weights  of  the  other  eigenwaveforms  eorresponding  to  other 
hypotheses  set  are  automatically  updated.  In  other  words,  the  excess  weight  from  the 
hypothesis  with  the  largest  (last)  likelihood  is  re-distributed  to  the  other  hypotheses. 

The  key  idea  in  improving  PWE  waveform  in  low  SNR  is  to  use  the  previous  and 
latest  weights  of  the  hypothesis  with  the  largest  current  likelihood  update  in  forming  the 
next  waveform  weights.  The  result  is  the  tempering  of  the  updates  whieh  swing 
dynamically  due  to  large  noise  realizations.  The  performanee  gain  of  this  teehnique  in 
low  SNR  over  MAP-PWE  and  MF-PWE  are  significant.  The  flow  diagram  of  the  TS- 
PWE  procedure  is  summarized  in  Fig.  30. 

Again,  we  set  up  a  Monte  Carlo  target  recognition  experiment  using  various  PWE 
waveforms.  In  Fig.  31,  mean  number  of  illuminations  performance  vs  transmit  energy  is 
shown  for  various  waveforms.  It  is  interesting  to  note  that  the  TS-PWE  is  better  than 
other  waveforms  in  low  transmit  energy  (less  than  -8  dB).  In  Fig.  32,  we  show  the 
identifieation  performance  results  of  the  MC  experiment  as  a  function  of  transmit  energy 
while  varying  the  number  of  transmissions  (or  updates).  Sinee  MF-PWE  and  TS-PWE  in 
Fig.  31,  we  utilize  these  two  in  Fig.  32.  Notiee  that  TS-PWE  outperforms  MF-PWE. 
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Flow  diagram/ procedural  steps  for  the  TS-PWE 


Figure  30.  Operational  flow  diagram  for  the  TS-PWE  scheme  with  a  closed-loop  radar. 
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Probability  of  Identification  p  Mean  nunnber  of  illuminations 


3 1 .  Mean  number  of  illuminations  as  a  function  of  (received)  transmit  energy  for 
various  adaptive  waveforms  assuming  unity  noise  energy. 


Ex(dB) 


Figure  32.  TS-PWE  vs  MF-PWE  identification  performance  comparison  as 
parameterized  by  the  number  of  updates  (or  number  of  transmissions)  assuming 

unity  noise  energy. 
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In  summary,  the  region  of  low  SNR  is  usually  the  most  interesting  in  any 
deteetion  problem  whether  it  is  radar,  eommunieations,  ete.  This  is  beeause  the 
probability  of  deteetion  is  a  funetion  of  SNR.  The  lower  the  SNR  means  the  lower  the 
probability  of  deteetion  where  there  is  usually  very  little  ean  be  done  exeept  the  obvious 
whieh  is  to  either  inerease  signal  energy  or  lower  noise  power.  Here,  although  both  MF- 
PWE  and  MAP-PWE  are  superior  to  various  waveforms  in  improving  probability  of 
deteetion  at  low  SNR,  we  are  still  able  to  produee  more  gain  with  TS-PWE  at  low  SNR. 
The  use  of  both  previous-to-last  and  latest  likelihoods  allows  for  the  TS-PWE  seheme  to 
temper  volatile  updates  due  to  large  noise  realizations.  As  sueh,  TS-PWE  performanee 
improvement  in  the  low  SNR  region  is  signifieant. 
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VIII.  MULTIPLE  TARGET  IDENTIFICATION  AND  RANGE- 

DOPPLER  LOCATION 


It  is  our  goal  to  simultaneously  identify  the  target  type,  determine  how  many 
targets  in  eaeh  type  and  eorreetly  loeate  the  ranges  and  Doppler  for  multiple  moving 
targets.  Range-Doppler  via  ambiguity  funetion  for  point  targets  is  well  eovered  in  the 
literature.  A  good  starting  point  is  [4].  In  this  seetion,  we  utilize  range-Doppler  map 
(RDM)  to  eorreetly  loeate  an  extended  target.  We  integrate  RDM  in  our  CR  platform  to 
simultaneously  identify  a  target  from  a  set  of  possibilities  with  the  use  of  our  adaptive 
PWE  waveforms.  We  investigate  different  scenarios  and  finally  form  our  comprehensive 
scheme. 

A,  TWO  TARGETS  OF  SAME  TYPE  SCENARIO 

In  this  section,  the  goal  is  to  jointly  locate  (in  range  and  Doppler)  and  identify  two 
targets  of  same  type  with  arbitrary  speeds  and  delays  (range).  Other  complicated 
scenarios  will  ensue  in  later  sections.  In  Fig.  33  we  illustrate  the  two  present  extended 
targets  in  RDM  where  a  small  amount  of  noise  is  added.  The  top  panel  is  a  three- 
dimensional  RDM.  The  bottom  panel  is  a  two-dimensional  RDM  (which  is  the 
conventional  way  of  illustrating  RDM).  For  simplicity  and  brevity,  we  choose  sample 
delay  (for  range)  and  normalized  frequency  for  Doppler  (for  velocity).  Both  parameters 
can  easily  be  converted  to  actual  range  and  velocity  from  the  RDM.  The  overall 
probability  of  correctly  locating  an  extended  target  in  range-and-Doppler  and  identifying 
target  type  is  clearly  a  function  of  received  signal-to-noise  ratio,  number  of  transmissions 
L ,  PWF  schemes  and  the  target  responses  themselves  via  maximum  target  eigenvalues. 
In  this  scenario,  we  assume  that  two  present  targets  are  of  the  same  type.  Since  we 
already  assumed  that  there  are  two  targets  present,  our  problem  therefore  is  to  figure  out 
where  these  two  targets  are  in  range-and-Doppler  and  identify  the  target  type  from 
various  possibilities  simultaneously.  In  practice,  a  priori  assumption  may  come  from 
other  sensors  or  some  intelligence  information.  To  perform  integrated  range-and-Doppler 
location  and  target  type  identification  given  they  are  of  same  type,  we  utilize  our  closed- 

loop  update  procedure  with  the  use  of  PWF  waveforms.  With  the  previous  identification 
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problem  where  transmission  number  is  fixed,  the  target  with  the  largest  weight  update  is 
chosen  to  be  the  detected  hypothesis.  Since  we  know  that  there  are  two  targets,  we 
choose  the  hypothesis  with  the  largest  weight  updates  after  a  fixed  number  of 
transmissions. 

We  now  describe  the  scenario.  Assuming  M  -  4  possible  extended  target  types, 
we  form  the  normalized  eigenwaveform  for  each  target  hypothesis  and  scale  each 
eigenwaveform  with  the  square-root  of  the  initial  waveform  weight  assigned  to  each 
target  and  then  form  the  first  PWE-based  waveform  (as  previously  when  we  assumed  the 
target  to  be  static).  The  difference  here  is  that  we  now  send  a  series  of  pulses  R  times  (in 
this  work,  R  =31),  since  we  are  interested  in  forming  RDMs.  For  the  received  return  from 
every  set  of  R  pulses  (recall  that  there  are  L  transmissions,  where  there  are  R  pulses  in 
each  transmission),  M  matched  filters  from  all  possible  extended  targets  are  applied  to 
form  the  Doppler  filter  banks  to  form  M  RDMs.  For  fair  comparison,  we  set  unity 
energy  of  each  transmission  (i?  pulses).  Assuming  N  is  the  length  of  target  impulse 
response,  the  matched  filtered  sequence  takes  on  a  length  of  4N  -  3  for  any  target 
present.  If  the  SNR  is  sufficient,  the  two  highest  magnitudes  in  the  RDM  may  indicate 
the  range -Doppler  locations  of  the  two  targets.  Due  to  noise  realization,  one  magnitude  is 
slightly  higher  than  the  other.  We  pick  the  4A-3  sequence  corresponding  to  the  larger 
one  to  use  for  our  likelihood  update  calculations.  To  properly  scale  the  next  PWF 
waveform  weights,  we  normalize  the  sum  of  the  likelihood  to  one  as  usual  and  form  the 
waveform  weights  such  that  the  energy  constraint  is  accommodated  (58).  We  continue 
until  we  send  the  fixed  number  of  transmission  L .  We  choose  the  target  type  with  the 
largest  likelihood  value  which  allows  us  to  pick  the  RDM  out  of  the  four  possible  RDMs. 
Since  there  are  two  targets,  we  pick  the  two  range-Doppler  cells  with  the  largest 
magnitudes.  Thus,  we  have  jointly  decided  the  target  type  and  located  the  range-Doppler 
cells  where  the  two  targets  may  be  located.  The  illustration  of  RDMs  updates  for  each 
target  type  is  shown  in  Fig.  33  to  Fig.  36  where  the  two  targets  can  be  easily  determined 
in  the  “decided”  RDM  after  five  updates  in  Fig.  36.  Here,  we  use  MF-PWF  as  the 
designated  PWF  scheme. 
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Figure  33.  Range  Doppler  map  of  two  targets  of  the  same  type  (some  noise  added). 


No  update  of  RDM  for  target  1 


No  update  of  RDM  for  target  2 


•ano* 

No  update  of  RDM  for  target  3 


1400 


No  update  of  RDM  for  target  4 


Figure  34.  Four  eandidate  RDMs  for  four  possible  target  types  (where  two  targets  of  the 
same  type  are  present)  with  PWE  after  one  transmission  (no  update). 
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Two  updates  of  RDM  for  target  1 


Two  updates  of  RDM  for  target  2 


Figure  35.  Four  candidate  RDMs  for  four  possible  target  types  (where  two  targets  of  the 
same  type  are  present)  with  PWE  after  three  transmissions  (two  updates). 
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Five  updates  of  RDM  for  target  3 
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Figure  36.  Four  candidate  RDMs  for  four  possible  target  types  (where  two  targets  of  the 
same  type  are  present)  with  PWE  after  six  transmissions  (five  updates). 
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In  Fig.  37,  the  probability  of  identification  {P.  which  is  the  probability  of 
correctly  determining  target  type),  which  is  the  probability  of  correctly  determining 
the  range-Doppler  location,  and  the  overall  or  joint  probability  (P)  with  one 

transmission  (no  update)  and  six  transmissions  (live  updates)  using  MF-PWE  scheme  are 
shown  via  MC  experiments.  The  overall  probability  is  the  probability  that  the  range- 
Doppler  location  and  target  type  are  both  correct.  In  other  words,  if  we  correctly 
determined  the  target  type  and  not  correctly  determined  the  correct  range-Doppler  cell 
(despite  the  fact  that  the  actual  RDM  peak  is  in  an  adjacent  cell),  we  still  count  this 
against  the  overall  probability.  Notice  also  that  the  probability  of  correct  identification  is 
better  than  the  probability  of  correctly  determining  the  range-Doppler  locations.  This  is 
because  in  our  experiment,  there  are  only  four  possible  target  types  while  there  are 
numerous  range-Doppler  cells.  We  can  improve  P^  by  increasing  the  range-Doppler  area 

(number  of  cells)  in  which  to  decide  where  the  target  may  be  located.  For  the  moment, 
this  is  outside  the  scope  of  our  study  but  definitely  is  an  interesting  one  to  tackle  at  a  later 
time. 

It  is  also  clear  from  Fig.  37  that  the  overall  performance  (bottom  panel)  is  jointly 
affected  by  probability  of  identification  and  probability  of  correctly  determining  the 
range  and  Doppler  location  (upper  and  middle  panels).  From  transmit  energy  of  -4  to  0 
dB  energy  units  which  lead  to  almost  perfect  identification  percentage  {P.=\),  P^  is 

almost  exactly  the  same  as  P^  .  Recall  that,  we  put  a  very  strict  definition  on  P^ ,  i.e.,  we 

only  declare  range-Doppler  detection  when  the  matched  filtered  peak  is  truly  located 
where  the  true  target  peak  would  be  (although,  the  peak  may  simply  be  in  the  adjacent 
range-Doppler)  cell.  It  is  also  important  to  notice  that  more  transmissions  help  both 
identification  rate  as  well  as  overall  performances. 
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Probability  of  correctly  determining  target  type 


Figure  37.  Probability  of  correctly  determining  target  type,  probability  of  correctly 

determining  the  range-Doppler  location,  and  the  overall  performanee  of  eorrectly 
determining  target  type  and  range-Doppler  loeation  using  MF-PWE  with  two 

targets  of  same  type. 


In  Fig.  38,  we  show  the  eomparison  of  identification  probability  (/’)  of  six 
updates  and  ten  updates  with  various  PWE  schemes  for  the  two  moving  extended  targets 
(of  same  type  scenario).  It  is  obvious  that  P.  is  function  of  both  PWE  method  and 

number  of  transmissions  (or  update  number).  It  is  clear  that  TS-PWE  has  the  best 
performance  in  terms  of  probability  of  correct  identification  especially  in  the  low-SNR 
area.  Since  ME-PWE  and  TS-PWE  are  again  the  best  performers,  we  utilize  both  in  Fig. 
39.  In  Fig.  39,  the  overall  probability  comparison  between  MF-PWE  and  TS-PWE  are 
shown  for  three  transmissions  and  six  transmissions.  Notice  that  for  three  transmissions, 
TS-PWE  performs  slightly  lower  than  ME-PWE.  However,  at  six  transmissions,  TS-PWE 
is  better.  This  means  that  although  TS-PWE  is  better  than  ME-PWE  at  target 
identification,  it  is  not  always  better  than  ME-PWE  in  probability  of  correctly 
determining  range-Doppler  location.  Higher  transmission  number  is  needed  for  that  to 
happen. 
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Probability  of  identification  (P)  of  various  PWE  schemes  with  multiple  updates 


Figure  38.  Identification  probability  comparison  of  six  updates  and  ten  updates  of  various 
PWE  schemes  assuming  unity  noise  energy  for  two  targets  of  same  type  scenario. 


Probability  of  correct  identification  and  location  with  PWE  schemes 


Figure  39.  Probability  comparison  of  correctly  determining  targets  type  and  range- 
Doppler  location  of  no  update,  one  update,  two  updates,  and  five  updates  vs  MF- 
PWE  and  TS-PWE  schemes  assuming  unity  noise  energy  with  two  targets  of 

same  type. 
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In  summary,  we  slightly  modified  the  CR  platform  to  adapt  the  assumption  of  two 
target  of  same  type  and  form  a  platform  that  ean  eorreetly  identify  the  target  and  eorreetly 
indicate  the  range-Doppler  location  or  cell  for  moving  extended  targets.  We  used  the 
three  PWE-based  waveforms  ah  of  which  are  effective  in  terms  of  three  probabilities 
mentioned  above.  TS-PWE  is  the  most  effective  in  terms  of  target  recognition 
specifically  in  low  SNR. 

B.  UNKNOWN  NUMBER  OF  TARGETS  (OF  THE  SAME  TYPE) 

In  this  section,  we  do  not  assume  the  number  of  targets  present  but  since  this  is  a 
target  recognition  problem  we  will  generate  at  least  one  target  present  in  our  MC 
experiments.  If  there  are  multiple  targets,  we  assume  they  are  of  the  same  type  in  this 
section.  Other  complicated  scenarios  are  tackled  in  later  sections.  We  will  utilize  ME- 
PWE  and  TS-PWE  schemes.  Here,  we  modify  the  receiver  to  incorporate  thresholds 
dictated  by  a  given  false  alarm  since  a  true  detection  problem  is  added  into  the  mix,  i.e., 
the  number  of  targets  present  is  not  known  (but  we  at  least  assume  one).  As  before,  the 
hypothesis  with  the  largest  updated  waveform  weight  is  deemed  to  be  the  true  target 
hypothesis.  The  sequences  containing  highest  magnitudes  in  the  chosen  RDM  may  be 
retrieved  for  target  range  and  Doppler  location.  The  question  here  is  to  determine  how 
many  targets  there  are  in  latest  updated  RDM  and  what  is  the  proper  threshold  to  be  used 
in  locating  targets’  range  and  Doppler. 

We  assume  targets  are  of  length  N  (in  this  work,  A  =  3 1 )  and  the  length  of  test 
sequence  y  containing  peak  value  in  RDM  is  4N  -  3  .  In  order  to  form  each  range 
Doppler  map,  R  pulses  are  transmitted  (in  this  work,  R  =  3 1 ).  The  H .  is  the  2  N  -by-  N 

convolution  matrix  of  target  type  i  and  w  is  the  complex -valued  AWGN  vector.  Each 
time  the  radar  transmits  a  waveform  x ,  a  noise-corrupted  version  of  the  reflected  target 
echo  is  received.  Since  we  have  four  possible  (i.e.,  different)  targets,  we  use  the  four 
different  matched  filters  as  usual  for  ME-PWE.  The  /y|,  (y)  which  is  the  likelihood  of 

hypothesis  given  target  is  of  this  form 

fj\i  (y)  =  P  exp[2  Re(sJ'y)]  (76) 
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where  we  ean  substitute  y  =  H,.x  + w  and  s^.  =  H^.x  into  (76). 


Before  we  ean  determine  range -Doppler  loeation,  we  have  to  eonsider  the 
problem  of  deteetion  for  the  entire  range-Doppler  map  since  we  do  not  know  the  number 
of  targets.  The  detection  hypotheses  are 


^0 :  y = w 
:y  =  H.x  +  w. 


(77) 


It  can  be  shown  that  the  natural  log  likelihood  for  the  null  hypotheses  given 
target  is  of  normal  distribution  of  the  form 


ln[/o|,  (y)]  =  In  /?  +  2  Re(sf  w), 


(78) 


and  the  hypotheses  is 

ln[yi|,  (y)]  =  In  +  2  Re(x"Hf  H,x,  +  x^Hf  w).  (79) 

It  can  be  shown  that  the  likelihood  ratio  threshold  r  for  given  probability  of  false 
alarm  and  AWGN  of  sample  variance  cr^  is  of  the  form 


'■=^2E,alQ-'(P,,)  +  lnp.  (80) 

Thus,  the  probability  of  a  target  being  in  a  specific  range  and  Doppler  location  for  target 
type  i  is 


Pn=Q{Q\PpA)  +  ,^} 

V  cr,. 


12^ 

2 


=  Q{Q\Ppa)+. 


2Tx"HfH,x 


=  Q{Q\Ppa)+ 


M 


M 


j=\  j=\ 


cr„ 


(81) 


where  the  transmission  number  L  is  incorporated  and  R.  is  the  autocorrelation  matrix  of 
target  response. 
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It  is  clear  from  a  previous  ehapter  that  is  not  unique.  The  value  of  is 

affeeted  by  transmit  signal  energy  and  weights  for  eaeh  eigenveetor,  target  response 
energy,  and  maximum  eigenvalue.  If  the  PWE  update  proeess  eonverges  to  the  true  target 
type  and  mostly  distributes  the  weight  to  target  type  i  (true  target),  the  deteetion 
probability  for  eaeh  range-Doppler  eell  given  eorreet  target  type  identifieation  is  shown 
in  (81).  However,  if  the  proeedure  makes  the  wrong  target  type  identifieation  (whieh 
means  most  energy  is  not  in  the  target  type),  the  probability  of  loeating  the  eorreet 
range-Doppler  eell  will  tend  to  be  very  low  as  shown  in  Fig.  40  where  target  1  is  the 
eorreet  target  type.  Wrong  identifieation  also  leads  to  poor  probability  of  loeating  range- 
Doppler  loeation. 

It  ean  be  shown  that  when  eorreet  target  type  identifieation  is  made  and  all 
transmit  energy  is  distributed  to  eorreet  target  type,  the  equation  (81)  ean  be  modified  as 


Po=Q{Q^\Pfa)+ 


/2F(q,"R,q,), 
_2  ^ 


<J.. 


=  Q{Q^\p,a)  +  , 


cr,! 


where  is  the  energy  of  target  response  and  is  assumed  to  be  unit  energy  in  this  seetion. 
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Probability  of  locating  correct  range  and  Doppler  31  Probability  of  locating  correct  range  and  Doppler 


;  40.  Probability  of  locating  correct  range  and  Doppler  given  target  type 
identification  and  probability  of  false  alarm  ( Pp^  =0.1)  assuming  unity  noise 
energy  where  type  1  is  the  correet  hypothesis. 


Probability  of  locating  correct  range  and  Doppler 
given  correct/incorrect  target  type  identification  and  various  probability  of  false  alarm 

_ 1  I  r 


.  PjB  =  0.1  (correct  idertificebon)  ' 
•  =  0.05  (correct  identificalion) 

-  P,e  =  0.001  (correct  identificationi 
.  P,a  =  0.1  (incorrect  idenbiication) 

■  P^  =  0.06  (incorrect  idenbbcabon) 

..  Pa  =  0.001  (incorrect  idenbbcabon)  i 


Ex(dB) 


Figure  4 1 .  Probability  of  locating  correct  range  and  Doppler  given  correct/incorrect 
target  type  identification  and  various  probability  of  false  alarm  assuming  unity 

noise  energy. 


79 


Also,  in  Fig.  41,  it  shows  clearly  that  the  probability  of  correctly  locating  range 
and  Doppler  is  intimately  tied  to  the  probability  of  target  type  identification.  The  plots  in 
Fig.  41  re-enforces  the  fact  that  once  a  wrong  identifieation  is  made  then  the  probability 
of  correctly  locating  the  range-Doppler  cell  of  a  target  becomes  mueh  lower. 

Again  we  set  up  Monte  Carlo  experiments  where  we  randomly  generate  the 
number  of  targets  and  the  range-Doppler  loeations  of  those  targets.  In  this  particular 
experiment  we  generate  one,  two,  or  three  targets  (of  the  same  type)  while  assigning 
range-Doppler  locations  (randomly)  in  each  experiment.  In  our  joint  recognition- 
detection  procedure,  reeall  that  we  identify  the  targets  by  ehoosing  the  RDM.  Again  since 
we  do  not  know  the  number  of  targets,  we  have  to  perform  detection  via  threshold  in  the 
entire  RDM.  The  threshold  is  based  on  the  PFA  given  a  target  type  from  (80).  Onee 
targets  are  “deteeted,”  then  we  try  to  determine  the  targets’  range-Doppler  loeations. 
Since  this  is  a  deteetion  problem,  the  threshold  used  dietates  the  aetual  .  Using  MF- 
PWE,  we  calculate  the  “strict”  overall  probability  P  (which  is  correct  identification, 
correct  number  of  targets  and  correct  range-Doppler  locations  all  at  the  same  time)  as  a 
function  of  transmit  energy.  Notice  that  is  tied  to  PFA  (since  P^  is  tied  to  PFA).  It  is 

clear  that  in  Fig.  42  that  the  threshold  (ealculated  from  probability  of  false  alarm)  and  the 
number  of  updates  affeet  overall  probability. 
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Joint  probability  of  correct  target  type  identification  and  range-Doppler  location  with  3  updates 


Figure  42.  Joint  probability  of  locating  correct  ranges,  Doppler  and  target  types  vs 
various  MF-PWE  updates  numbers  given  probability  of  false  alarm  assuming 

unity  noise  energy. 


C.  TWO  TARGETS  OF  DIFFERENT  TYPES  (SCENARIO) 

In  this  section,  we  make  the  problem  a  little  more  interesting.  Here  we  allow  for 
different  target  types  to  be  in  the  scene.  For  now,  let  us  limit  the  number  of  targets  to  two 
(and  then  weTl  consider  more  in  a  later  section).  Assuming  the  number  of  targets  (here 
two)  removes  the  detection  component  of  the  signal  processing  since  we  already  know 
the  number  of  targets.  This  is  akin  to  the  “two  targets  of  the  same  type  scenario.” 
However,  in  this  problem,  we  have  two  different  targets. 

We  may  consider  the  problem  that  of  a  multiple  hypothesis  testing  (MHT) 
problem.  Unfortunately,  in  practice  when  there  are  more  possible  targets  to  consider  and 
when  there  are  more  scenarios  to  consider  (e.g.,  ten  targets  total,  two  targets  are  type  1,  2 
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targets  are  type  2,  one  target  is  type  3,  etc.),  then  the  number  of  hypotheses  increase 
dramatically  due  to  the  combinatorial  implications  to  the  number  of  hypotheses.  If  we  are 
to  field  a  practical  cognitive  radar  system,  then  its  computing  resource  has  to  be 
reasonable.  As  such,  we  will  focus  on  how  to  update  the  waveform  weights  via  the 
likelihoods  corresponding  to  the  original  four  hypotheses  in  the  target  recognition 
problem.  This  reasonably  makes  sense.  For  example,  if  there  are  two  different  target 
types  (say  type  1  and  type  3)  in  four,  it  would  make  sense  that  after  a  few  transmissions 
that  the  Wj  and  approach  0.5  such  that  the  overall  waveform  favors  those  two  types. 

We  illustrate  the  procedure  in  Fig.  43.  Amazingly,  there’s  very  little  difference 
between  this  procedure  to  that  of  “two  targets  of  the  same  type  scenario.”  Here, 
effectively  the  only  difference  is  that  we  choose  two  RDMs  in  the  end.  Also,  the 
waveform  weights  will  be  distributed  to  the  two  target  types  (about  0.5  each  when 
homing  in  on  the  two  targets). 

We  conduct  an  experiment  where  two  targets  are  present  (type  1  and  type  4)  using 
the  procedure  in  Fig.  43.  We  show  the  4  RDMs  after  the  first  transmission  (no  waveform 
update)  in  Fig.  44.  Notice  that  it  is  difficult  to  tell  where  the  targets  are.  Now  we  show 
what  happens  to  the  RDM  after  seven  transmissions  in  Fig.  45.  Notice  that  the  two  peaks 
in  the  two  RDMs  corresponding  to  targets  type  1  and  type  4.  Clearly,  the  procedure  can 
correctly  identify  two  targets  and  correctly  determine  range -Doppler  locations. 
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Flow  diagram /procedural  steps  for  the  two  different  targets 


Decide  two  target  types  and  locate  range-Doppler  cells 


Figure  43.  Flow  diagram/procedural  steps  for  two  different  target  types  (seenario). 
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Range  Doppler  Map  after  matched  filter  1 


Range  Doppler  Map  after  matched  filter  2 


SMI 


Ring*  (SiHnpl») 


<00 


Range  Doppler  Map  after  matched  filter  3 


Range  Doppler  Map  after  matched  filter  4 


80- 

HI- 


Rkng*  (SatnplA) 


Figure  44.  No  update:  RDMs  of  two  targets  from  different  types  after  M  matched  filters 

(M  =  4). 


Range  Doppler  Map  after  matched  filter  1 
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Range  Doppler  Map  after  matched  filter  3 
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Range  Doppler  Map  after  matched  filter  2 
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Range  Doppler  Map  after  matched  filter  4 
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Figure  45.  Six  updates:  RDMs  of  two  targets  from  different  types  after  M  matched  filters 

(M  =  4). 
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Via  Monte  Carlo  experiments,  we  illustrate  the  overall  (or  joint)  probability  of 
correctly  identifying  (both  target  present)  and  their  range-Doppler  locations  in  Fig.  46 
with  the  use  of  MF-PWE  and  TS-PWE  adaptive  waveforms  as  a  function  of  transmit 
energy  while  varying  the  number  or  transmissions  (or  updates).  If  we  compare  Eig.  39 
(same  target  type)  and  Eig.  46  (different  target  types),  it  is  clear  that  performances  are 
lower  for  Eig.  46.  This  is  because  the  waveform  weights  are  distributed  between  two 
targets  (for  different  types)  while  the  weight  gets  distributed  mostly  to  one  target  (for 
same  target  type),  which  seems  intuitive. 


n 

(Q 

0 

CL 


Ex(dB) 


Eigure  46.  Overall  probability  of  identification  and  location  for  two  targets  of  different 

types  assuming  unity  noise  energy. 


D.  UNKNOWN  NUMBER  OF  TARGETS  AND  DIFFERENT  TYPES  OF 
TARGETS 

In  this  section,  we  assume  that  the  number  of  targets  is  not  known  or  there  could 
be  multiple  targets  in  the  scenario.  Moreover,  there  could  be  multiple  target  types. 
However,  we  assume  that  number  of  target  types  present  is  less  than  the  possible  target 
types.  This  assumption  is  from  the  intuition  gained  from  the  previous  section.  If  there  are 


85 


two  out  of  four  target  types  then  the  waveform  weights  for  the  two  targets  will  eventually 
reaeh  -0.5.  If  there  are  three  out  of  four,  then  the  weights  will  eventually  reaeh  -0.33.  So 
for  four  out  of  four,  then  it  will  simply  revert  to  0.25  as  intuition  would  dietate.  Our  goal 
is  to  utilize  the  insights  gained  from  the  three  scenarios  above  to  form  a  comprehensive 
scheme  to  simultaneously  figure  out  the  number  of  targets,  identify  target  types,  and 
determine  the  range-Doppler  cells  of  these  targets. 

Let  us  start  from  a  general  scenario  where  there  are  M  possible  target  types 
(regardless  of  the  number  of  targets  present).  If  there  is  one  target  type,  then  the  large  or 
“high”  waveform  weight  distribution  goes  to  the  correct  target  hypothesis  (assuming 
enough  SNR  and  enough  transmissions)  which  approaches  1  as  the  platform  homes  in  on 
the  correct  target  type  while  the  other  (“low”)  weights  approach  zero.  For  two  target 
types,  the  “high”  weights  would  approach  0.5  (assuming  enough  SNR  and  enough 
transmissions)  and  so  on  as  so  forth.  As  the  number  M  of  the  possible  types  increases, 
the  weight  differences  between  the  group  of  “high”  weights  and  the  group  of  “low” 
weight  decrease.  In  other  words,  it  may  be  harder  to  correctly  identify  how  many  types  of 
present  targets  there  are  as  the  number  of  present  target  types  increases.  It  can  be  shown 
that  the  weight  distribution  for  MF-PWE  for  M  possible  target  types  where  j  indicates 
targets  that  are  not  present  and  i  indicates  targets  that  are  present  is  the  following; 


W,= 


EU,,{x)-\ 

V/  y/ 


(82) 


where  ^  1’  2, .., i  index  for  present  target  type.  In  (82),  recall  that  we  assume  (x 

for  simplicity. 

The  first  term  in  the  denominator  is 
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(83) 


r  M  M 

Z  £[/;„  (j:  )] = Z  /» (Z  >"  R, 

Vi  Vi  -^1  r=l  r=l 

which  is  the  summation  of  all  likelihoods  of  present  targets.  And  the  seeond  term  is  the 
summation  of  all  likelihoods  of  targets  not  present  and  is  given  by 

'•y  P  M  M 

Z^[/>i.(*>l  =  Z/’«>‘P[-^(ZVkq.>"H“H,(ZVii;:q.)].  (84) 

vy  Vi  ^1  r=\  r=\ 

The  TS-PWE  weight  distribution  of  M  possible  types  is  more  diffieult  to  analyze 
beeause  it  involves  the  modified  likelihood  funetions. 

In  this  final  seetion,  we  also  inelude  the  ease  where  there  may  not  be  a  target.  In 
seetion  B,  although  we  assumed  the  number  of  targets  to  be  unknown,  we  did  assume  in 
our  experiments  that  at  least  one  target  is  present.  In  other  words,  the  set  of  algorithms  in 
seetion  B  has  to  be  modified  to  aeeommodate  the  general  case  when  there  is  no  target 
present.  So  here  our  strategy  is  to  first  figure  out  if  there  is  a  target  or  not  in  the  seenario 
(after  a  number  of  transmissions/updates  say  L  ).  If  there  is,  then  we  figure  what  target 
type  with  another  set  of  transmission  (same  as  before  i.e.,  L),  and  then  another  until  our 
algorithm  says  there  is  no  more  target  (type)  present.  To  this  end,  the  target  identifieation 
proeedure  is  modified  by  adding  a  user-defined  threshold  if/  (a  pereentage)  and 
eomparing  that  to  a  formulated  measured  y/  from  the  latest  weight  distribution.  The 
threshold  iff  as  in  any  threshold  maybe  adjusted  to  improve  deeision-making  of  the 
algorithm.  The  smaller  the  threshold  iff  is  defined,  the  more  sensitive  or  the  more 
probable  it  is  to  identify  and  deteet  present  targets. 

If  there  are  M  possible  targets,  then  the  “high”  weights  that  will  be  approaehed 
by  our  likelihood  update  procedure  is  100  /  M  .  The  measured  iff  is  then 

iff  =  E{wJ- (85) 

where  are  the  weights  that  are  above  the  (100 /M)  and  are  lower.  If  iff  is  less 
than  the  desired  iff ,  then  the  threshold  is  not  erossed  and  thus  a  target  type  is  not 
deteeted.  If  iff  is  greater  than  designated  iff ,  then  we  look  for  the  target  type  that  has  the 
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highest  weight  and  deem  that  to  be  a  present  target  type.  Sinee  it  is  deemed  present,  the 
eigenwaveform  eorresponding  to  that  target  is  now  removed  from  the  next  L  number  of 
transmissions.  Now  we  normalize  the  remaining  weights  as  dictated  by  (57,  61)  and  try  to 
ascertain  if  there  are  remaining  targets.  In  other  words,  we  terminate  the  transmissions 
when  y/  is  not  greater  than  our  threshold  y/  . 

The  full  procedure  for  an  experiment  is  summarized  in  Fig.  48  and  enumerated 

here. 

1.  Assuming  M  possible  extended  targets,  generate  the  normalized 
eigenwaveform  for  each  target.  Decide  which  PWE  waveform  is  used  and  set  the  user- 
defined  threshold  y/ .  Scale  each  eigenwaveform  with  the  square-root  of  the  initial  weight 
assigned  to  each  hypothesis  and  then  form  the  first  PWE-based  waveform. 

2.  After  the  first  return,  calculate  the  likelihood  according  to  the  choice  of  PWE 
scheme.  Eor  MAP -PWE,  use  (59)  and  (69)  is  for  MF-PWE  and  TS-PWE.  Continue 
transmissions  until  the  desired  number  of  transmissions  L  (or  updates  which  is  L-l). 
Note  here,  the  TS-PWE  requires  at  least  2  transmissions  to  implement. 

3.  Update  the  weight  distribution  (75)  and  calculate  the  measured  y/  by  (85). 

4.  Decide  if  there  is  a  target.  Recall  that  target  present  is  decided  if  the  measured 
y/  is  greater  than  designated  y/ .  In  our  experiment,  y/  is  set  10  percent.  If  y/  is  less  than 
designated  value,  go  to  step  7. 

5.  Determine  the  number  of  targets  from  the  latest  RDM  of  the  decided  type.  This 
is  performed  by  using  the  traditional  threshold  to  the  entire  RDM.  The  threshold  is 
calculated  from  (80)  and  is  a  function  of  probability  of  false  alarm. 

6.  Go  back  to  step  1  and  eliminate  the  accepted  target  type.  That  is,  M  -  M  -I . 
Also,  we  remove  the  likelihood  for  this  target  type  in  calculating  the  new  waveform 
weights  as  dictated  by  (56,  57). 

7.  If  is  less  than  yr ,  decide  there  is  no  more  target  present  and  end  the 
algorithm. 
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The  overall  (or  joint)  performanee  of  eorreetly  identifying  the  number  of  targets, 
target  types,  and  loeating  range-Doppler  eell  for  those  targets  is  shown  in  Fig.  47.  In  this 
work,  we  use  Af  =  4 ,  y/  =  \Q  and  P^^=0.1  to  eompare  with  various  waveforms.  Sinee  our 

proeedure  identifies  each  target  every  L  number  of  transmission  (or  L-\  updates),  the 
plots  are  parameterized  by  number  of  updates  per  target  type.  It  is  obvious  that  both  the 
MF-PWE  and  TS-PWE  perform  much  better  than  wideband  waveform  as  may  be 
expected.  The  algorithm  presented  here  (used  in  a  CR  platform)  clearly  can  be  used  for 
the  general  problem  of  figuring  out  the  number  of  target  present,  target  type 
identification,  and  correctly  locating  range-Doppler  cells  of  those  targets.  The  flow 
diagram  is  summarized  in  Fig.  48. 


Figure  47.  Overall  probability  of  identification  and  location  for  multiple  targets  of 
unknown  types  assuming  unity  noise  energy  (^^  =  10  percent  and  Pp^  =  0.1). 
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Figure  48.  Flow  diagram/procedural  steps  of  multiple  targets  identifieation  and  range- 

Doppler  loeation  algorithm. 
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IX.  CONCLUSION 


In  this  work,  we  set  out  to  fulfill  three  main  related  objeetives.  The  first  is  to 
investigate  the  applieation  of  eigenwaveform  to  ambiguity  function  with  the  intent  of 
integrating  PWE-based  waveforms  with  range-Doppler  techniques.  The  second  objective 
is  to  extend  PWE-based  waveforms  and  produce  other  waveforms  based  on  PWE  for 
identification  of  extended  targets.  And  finally,  the  third  objective  is  to  design  an 
integrated  scheme  for  the  combined  problem  of  range-Doppler  location/detection  with 
extended  target  type  identification  with  the  use  of  a  closed-loop  radar  system  employing 
adaptive  waveforms. 

We  successfully  investigated  eigenwaveform-based  ambiguity  functions  with  the 
use  of  three  types  of  targets  (oscillatory,  non-oscillatory,  and  one  based  of  a  Boeing 
aircraft).  We  found  that  the  AE  properties  (peak  and  volume)  for  eigenwaveforms  to  be 
much  larger  than  that  of  traditional  waveforms  for  point  targets.  Ambiguity  function  is 
the  basis  of  range-Doppler  map  technique.  Various  range-Doppler  trade-offs  were  made 
with  the  use  of  pulsed  eigenwaveforms. 

We  made  improvements  to  the  MAP-PWE  by  proposing  two  new  PWE-based 
waveforms;  ME-PWE  and  TS-PWE.  Both  waveforms  performed  well  in  terms  of 
classification  over  MAP-PWE  and  classical  wideband  waveform. 

Ultimately,  we  used  the  PWE-based  waveforms  and  range-Doppler  techniques  to 
form  an  integrated  target  identification  and  range-Doppler  target  location  for  moving 
extended  targets.  Our  integrated  scheme  with  the  use  of  PWE-based  waveforms  worked 
well  against  the  traditional  wideband  waveform  in  terms  of  joint  probability  of  detecting 
targets,  identifying  target  types,  and  determining  range-Doppler  locations  for  those 
targets. 
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APPENDIX 


In  this  appendix,  we  diseuss  how  we  implemented  our  own  ambiguity  funetion  (in 
MATLAB)  and  how  we  verified  our  results. 

Let  x[/]  be  a  eomplex -valued  diserete-time  signal  that  we  would  like  to  analyze 
with  the  ambiguity  funetion.  For  compaet  presentation,  x  is  the  veetor  notation  of  x[/] . 
N  is  the  length  of  veetor  x  whieh  means  x{i\  contains  values  only  when 
i  =  0,1,2...(A-1)  .  We  will  interchange  these  notations  in  the  equation. 


The  definition  of  discrete  ambiguity  function  is  given  by 


^x[i]xx*[i-TjxW, 


fdn^i 

N 


(=0 


(1) 


where  is  the  delay  index  and  /^„is  the  Doppler  spread.  Note  here  0,  1,  ...,  N -\ 
and 

jlx— 

=  e'  ^  . 

In  general  the  length  of  may  be  desired  to  be  different  from  A .  For  now,  we  keep 
them  the  same  and  will  expand  it  later. 

Now,  we  keep  the  Doppler  spread  index  fixed  in  (1)  and  vary  the  delay  index 
for  comparison.  For  compact  presentation,  only  the  j  function  is  listed  out.  The 
equations  for  positive  are 
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;r[2;  /,„]  =  X 

i=2 

M3;/rfJ  =  ^x(zV(z 

i=2 
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li  x(zVO'-A^+iy^"^^”"'. 
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And  the  equations  for  negative  are 


N-2 


Ziyfdn]  = 

/=o 

Zi-2;  /dJ  =  Z  +  2)6^"'^^"""' 

z=0 
A/ -4 
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M-(iV-l);/,J  =  (/  + 

i=0 

Interestingly,  there  is  a  systematic  rule  for  the  j  function  for  fixed  Doppler 
spread  index  .  That  is 


M-(iV -!);/.„] 


M-2;/.J 

M0;/,J 
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where  the  K  matrix  is  defined  as 


x*(A-l) 

0  0 
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0 
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0 

0  0 

0 

x(0) 

which  is  taking  the  reversed  and  conjugated  x  and  is  lined  up  in  a  special  format  similar 
to  the  convolution  matrix.  And  the  f  matrix  is  defined  as 


x(2)  X 


(6) 


This  is  mathematical  vector  for  fixed  Doppler  spread  . 

For  general  ambiguity  function  in  matrix  form  (which  means  all  discrete  Doppler 
spreads  are  considered),  it  is  simply  the  expansion  of  (4)  as  in  following: 


95 


z[-iN-V);fd(M 

M-2;/.o] 

... 

Z\.~'^:'  fd{M-l)\ 

... 

Z\.~\  fd(M-\)\ 

MO;/,o] 

z[0;fd,] 

Z\S^->fd(M-l)^ 

ZW->fd(M-\)^ 

z[2;fdo] 

z[2;fdA 

/f[2?y^(M-i)] 

z[N-l;fdo] 

z[N-i;fd,] 

Z\.N  ~Xfd{M-\ 

(V) 


(2N-l)xN 


[F] 


NxM 


where  we  expand  the  one  fixed  Doppler  spread  to  all  running  indiees  by  modifying  f 
matrix  to  F  matrix  to  include  all  the  Doppler  spread  indices.  The  F  matrix  is  given  by 


[F], 
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(8) 


Therefore,  the  x  matrix  can  be  calculated  from  the  matrix  operation  of 


I  |2 

and  the  ambiguity  function  is  \x\j„\fdn\  ■ 

It  is  necessary  to  add  a  constraint  for  M  >  N  to  maintain  the  integrity  of 
ambiguity  matrix.  The  physical  interpretation  is  the  number  of  Doppler  components 
needs  to  be  at  least  as  much  as  the  length  of  x . 

96 


The  accuracy  of  operation  is  easy  to  check  from  definitions  of  zero-delay  cut  and 
zero-Doppler  cut.  By  definition,  the  zero-delay  cut  is  exactly  the  inverse  Fourier 
transform  of  x\i^  and  the  zero-Doppler  cut  is  the  autocorrelation  function  of  x\i\ .  It  can 
be  shown  from  the  figure  below  where  the  zero-delay  cut  and  zero-Doppler  cut  are 
exactly  the  same  as  the  autocorrelation  and  IFFT  results  using  Matlab’s  own  corr  and  fft 
functions. 


Zero-Doppler  cut 


Sample 


Zero-Delay  cut 


IFFT  result  from  x[n]^ 


Figure  1;  Comparison  of  ambiguity  function  cuts  to  theoretical  function. 
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